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I.  Identification  and  Significance  of  the  Problem 


The  nuclear  analog  of  the  ruby  lasep  would  embody  the  most 
straightforward  of  the  several  concepts  for  a  gamma  ray  laser. 
Not  surprisingly,  the  greatest  progress  toward  an  ultrashort 
wavelength  laser  has  been  achieved  along  that  direction  . 

For  ruby  the  identification  and  exploitation  of  a  bandwidth 

funnel  were  the  critical  keys  in  the  development  of  the  first 

laser.  There  was  a  broad  absorption  band  linked  through  efficient 

cascading  to  the  narrow  laser  level.  The  search  for  the  analogous 

situation  on  the  nuclear  scale  is  proceeding  quite  fruitfully. 

Significant  amounts  of  nuclear  fluorescence  have  been  pumped  by 

pulsed  bremsstrahlung  sources;  thus  proving  that  the  concept  of 

bandwidth  tunneling  works  as  well  at  the  nuclear  level  as  it 

does  for  the  ruby^.  This  critical  process  was  found  to  produce 

eleven  orders  of  magnitude  increase  in  the  amount  of  x-ray  pump 

77 

energy  absorbed  into  demonstration  nuclei  of  Se. 

Most  recently,  the  first  of  the  29  candidate  materials  for  a 
gamma  ray  laser  was  successfully  pumped  with  intense  pulses 
of  bremsstrahlung.  The  yield  of  nuclear  fluorescence  was  about 
10,000  times  greater  than  what  was  predicted  by  even  the  theories 
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including  the  large  factor  for  funneling  .  Selected  over  28 
alternatives  solely  on  the  basis  of  availability,  this  first 
material  narrowly  missed  the  threshold  of  feasibility  for  a  gamma 
ray  laser.  These  results  with  a  relatively  unattractive  candidate 
indicate  that  the  probabilities  should  be  raised  for  full  success 
of  one  of  the  other  28  materials. 

While  the  rate  of  progress  is  very  high  in  the  existing  SDIO 
Gamma  Ray  Laser  Program,  pacing  constraints  arise  from: 

1)  limitations  of  supporting  technology,  and 
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2)  extreme  scarcity  of  the  29  candidate  materials. 
Successes  in  alleviating  the  first  problem  would  make  it  possible 
to  proceed  more  rapidly  without  having  to  overcome  the  much  more 
expensive  aspects  of  the  second  problem.  Once  more  advanced 
technology  is  available  to  identify  the  best  of  the  29  from 
minute  samples,  the  problems  of  the  production  of  useful  amounts 
need  only  be  solved  for  the  one  material  selected.  It  was  the 
purpose  of  this  SBIR  effort  to  explore  the  development  of  one 
innovation  which  can  materially  advance  the  technology  supporting 
the  Gamma  Ray  Laser  Program. 

A  number  of  fundamental  factors  conspire  to  make  it 
extremely  difficult  to  accurately  measure  the  spectral  content  of 
intense  pulses  of  bremsstrahlung  containing  photon  energies  above 
200  keV.  However  this  is  just  the  type  of  pump  radiation 
currently  used  in  testing  materials  which  might  serve  as  nuclear 
analogs  of  the  ruby  system.  With  no  quantitative  measures  of 
spectral  input,  determinations  of  fluorescence  efficiencies  are 
extremely  difficult  and  even  the  relative  rankings  of  the  merits 
of  specific  nuclei  must  be  equivocal. 

There  are  no  periodic  spacings  of  atoms  or  lattices 
comparable  to  a  wavelength  at  such  energies  so  dispersion  is 
ineffective.  In  pulsed  sources  the  great  number  of  individual 
photons  insures  pile-up  in  solid  state  detectors  effectively 
eliminating  spectral  information.  The  attenuation  of  radiation  in 
matter  is  no  longer  structured  at  such  energies,  so  neither 
selective  filtering  nor  even  bandpass  layers  can  be  effectively 
configured.  Difficulties  are  so  severe  that  bremsstrahlung 
sources  of  intense  pulses  are  routinely  "calibrated"  by 
calculating  what  the  intensity  should  be.  Until  this  year,  this 
had  been  the  only  available  procedure  for  both  medical  linacs^ 
and  nuclear  simulators^. 
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A  very  recent  spinoff  of  the  SDIO  Gamma  Ray  Laser  Program  at 

the  University  of  Texas  at  Dallas  has  been  the  technique  of  X-ray 

7  8 

Activation  of  Nuclei  '  ( XAN )  for  the  calibration  of  intense 

bremsstrahlung  sources  operating  in  the  range  0.2  to  1.4  MeV.  As 

..  ,  .  ,  .  77  7  9  111 _  -  1  1  5  _ 

currently  implemented,  five  isotopes,  Se,  Br,  Cd,  In, 

and  1 ^7Er  are  used  to  sample  narrow  spectral  slices  of  the 

bremsstrahlung  fluence  illuminating  a  target.  The  isotopes  are 

excited  to  long  lived  states  whose  nuclear  fluorescences  can  be 

"read-out"  at  a  later  time  at  which  there  is  less  electrical 

noise.  With  five  isotopes  the  XAN  method  accommodates  measurement 

at  five  photon  energies  in  the  range  0.2  to  1.4  MeV. 

One  of  the  unexpected  results  of  the  pumping  of  the  first 

9 

actual  laser  candidate  was  the  importance  of  nuclear  states 

lying  above  1.4  MeV.  Unfortunately,  these  levels  could  only  be 

accessed  by  the  uncalibrated  parts  of  the  source  spectrum  for 

which  XAN  cannot  be  used.  The  technical  need  that  was  addressed 

by  the  research  of  this  Phase  I  effort  was  the  absolute 

1  0 

calibration  of  the  rest  of  the  range  of  photon  energies,  1.4  to 
6.0  MeV,  which  may  be  useful  in  studying  the  candidates.  In 
principle  this  can  be  done  with  an  extension  of  the  XAN 
instrumentation,  but  the  initial  calibration  and  validation  will 
require  technology  not  presently  available  to  the  existing  gamma 
ray  laser  projects  nor  currently  budgeted  for  development.  The 
research  we  have  completed  under  Phase  I  of  the  SBIR  program 
has  proven  that  the  integration  and  use  of  the  requisite 
technology  is  definitely  feasible. 


II.  Technical  Background 


Energetics  of  the  XAN  procedures  are  shown  in  Fig.  1.  As 
seen  there,  the  resonant  absorption  of  a  bremsstrahlung  photon, 
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Figure  1  :  Schematic  representation  of  the  energetics  of  the 
XAN  process  for  sampling  x-ray  spectral  intensities. 
Absorption  of  an  x-ray  first  excites  a  short-lived  state 
which  then  decays  in  part  by  cascading  to  a  long-lived 
state. 

symbolized  by  the  upward  arrow,  raises  a  nucleus  from  the  ground 
state  to  a  funneling  level  represented  by  the  partly  shaded 
rectangle.  Its  extremely  short  lifetime  gives  it  a  relatively 
broad  width  for  the  absorption  of  a  photon  from  the  pump 
continuum.  However,  broad  at  the  nuclear  level  is  still  narrow  in 
comparison  to  any  structure  in  the  x-rays  from  the  pump,  so  the 
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excitation  of  the  funneling  level  occurs  to  an  extent  dependent 
upon  the  intensity  in  a  narrow  slice  of  the  spectrum  of  the  input 
radiation.  Nonresonant  contributions  to  such  excitations  have 
been  shown  to  be  negligible  '  so  the  population  of  the 
funneling  level  will  be  simply  proportional  to  the  input 
intensity  at  the  corresponding  energy. 

As  shown  in  Fig.  1  by  the  shaded  portion  of  the  rectangle,  a 
fraction  bQ  of  the  population  cascades  to  the  longer-lived 
isomeric  state  which  ultimately  radiates  the  output  fluorescence. 
The  complementary  fraction  b  decays  back  to  the  ground  state. 
The  quantitative  yield  of  fluorescence  can  be  written  in  terms  of 
the  spectral  fluence  input  and  the  nuclear  parameters  of  Fig.  1 
in  a  particularly  convenient  form  here  reproduced  from  Ref.  13 
with  permission  of  the  authors. 


As  described  in  the  article  reporting  the  previous  implementation  of  this 
technique  at  lower  energies,  the  number  of  excited  nuclei  produced  in  the 
course  of  an  irradiation  of  the  type  shown  in  Fig.  1  can  be  computed  from  the 
expression: 


(*bab0oer/2),  *(Ej) 


(1) 


where  N  is  the  number  of  absorbing  nuclei  in  the  sample  and  the  summation  is 
taken  over  the  properties  and  pump  fluxes  characteristic  of  each  of  the  i-th 
broad  pump  bands  centered  at  pump  energy  E,-.  Only  one  such  pump  band  is  shown 
in  the  scheme  of  Fig.  1,  but  there  could  be  several  at  different  each 
funneling  its  population  into  the  same  output  level. 

In  Eq.  (1)  the  first  ratio  in  the  summation  is  composed  of  the  nuclear 
parameters,  while  the  second  describes  the  intensities  of  the  pump  x-rays 
which  are  assumed  to  be  continuous,  at  least  without  structure  on  the  fine 
scale  of  the  nuclear  absorption.  In  particular,  the  combination  ^(Ej)/A  is 
the  spectral  fluence  at  the  energy  E{  in  units  of  keV/keV/cm2.  Tacitly,  it 
has  been  assumed  that  the  duration  of  the  pump  source  is  less  than  the 
fluorescence  lifetime,  r0.  7 


Of  the  nuclear  parameters,  T  is  the  natural  width  in  keV  of  the  i-th  pump 


band,  as  shown  in  Fig.  1, 

r-*Ap 


(2) 


and  the  branching  ratios,  b#  and  b0,  give  the  probabilities  for  the  decay  of 
the  broad  level  back  into  the  initial  and  fluorescence  level,  respectively. 
The  pump  energy,  Es ,  is  in  keV  and  c0  is  the  Breit-Vigner  cross  section  for 
the  absorption  transition, 


A2  2Ie+l  1 
1*  2Ig+l  Qp+1 


(3a) 


where  A  is  the  wavelength  in  cm  of  the  gamma  ray  at  the  resonant  energy,  E, ; 
Ie  and  Ig  are  the  nuclear  spins  of  the  excited  and  ground  states,  respective¬ 
ly;  and  ap  is  the  total  internal  conversion  coefficient  for  the  broad  level  of 

Fig.  1. 


The  combination  appearing  in  parentheses  in  the  numerator  of  the  "nuclear 
part"  of  Eq.  (1)  is  the  integrated  cross  section  which  is  conveniently 
measured  in  units  of  10-29  cm2  kev, 


q,  -  (mbab0^or/2)i  .  (3b) 

This  expression  is  most  often  used  in  numerical  computation,  followed  closely 
by , 


Q,  -  qj/Ej  .  (3c) 

and  the  effective  partial  width  for  excitation  of  the  output  channel, 

-  (b,bor)f  .  (3d) 


On  this  scale  for  integrated  cross  sections,  unity 

or  even  a  few  tenths  suffice  for  a  material  to  be  useful  in 
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to  1 0 


12 


14 


sampling  spectral  f luences  of  the  order  of  10 

keV/keV/cm  at  the  corresponding  energy.  These  integrated  cross 

sections  are  the  critical  numbers  which  must  be  determined  before 

XAN  can  be  used  for  calibrating  bremsstrahlung  sources.  For 

2 

example,  the  entire  database  which  currently  supports  use  of  the 

XAN  process  for  calibrating  nuclear  simulators  over  the  0.2  to 

1.4  MeV  range  of  energies  is  shown  in  Table  I.  It  is  the 

extension  of  this  tabulation  to  a  higher  range  of  energies  that 
comprised  the  larger  objective  of  this  work.  How  to  accomplish 
this  goal  was  the  particular  subject  of  this  effort. 


1  4 

It  has  already  been  established  that  even  the  most  current 

nuclear  database1'’  of  general  use  does  not  offer  significant 

precision  for  nuclear  parameters  such  as  bQ  b^  to  allow  the  cross 

sections  to  be  calculated  to  a  useful  accuracy.  To  determine  them 

currently  requires  the  application  of  Eq.  (1)  to  sets  of 

fluorescence  data  actually  measured.  It  can  be  readily 

appreciated  that  sequences  of  measurements  of  the  numbers  of 

isomeric  nuclei,  Sm  produced  in  the  m-th  particular  irradiation 

delivering  fiuences  of  F  (E.)  at  each  of  the  energies  E. 

mi  l 

corresponding  to  a  gateway  cascading  to  the  isomer  lead  to  a  set 
of  equations, 


S1  =  Q1F1 (E1 )  +  Q2F1 (E2)  + 

s2  =  q1f2(e1)  +  q2f2(e2) 


,  (4a) 

,  (4b) 

,  (4c) 


where  the  are  the  combinations  of  nuclear  parameters 
comprising  the  first  ratio  in  Eq.  (1)  for  the  gateway  at  E^.  In 
principle,  the  technique  for  determining  the  set  of  consists 
of  measuring  the  Sm  obtained  with  a  number  of  different  spectral 
distributions  F  (E^)  and  then  solving  the  set  of  equations  for 
the  Ch  .  In  practice,  the  problem  is  much  more  difficult.  The 
gateway  energies,  E^  are  not  usually  known,  and  the  problem  of 
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Table  I 


Summary  of 

nuclear 

parameters  for 

applications  of  the 

Activation 

of  Nuclei 

at  energies  below 

1 . 4  MeV. 

Nuclide 

E  . 
m 

Eout 

Integrated 

Cross  Section 

(keV) 

(keV) 

-  2  9  2 

(10  cm  keV) 

79 

Br 

761 

207 

6.2 

77  „ 

Se 

250 

162 

0.20 

4803 

0.87 

81  8 

0.7 

1005 

30.0 

115- 

In 

1  078 

337 

20.0 

1  1 1  _  . 

Cd 

1  1  90 

245 

9.8 

1  67^ 

Er 

900 

208 

23.0 

aThe  effects  of  the  440  and  521  keV  transitions  have  been 
combined. 


determining  the  F^E.^)  is  the  general  one  motivating  the 
measurements  of  the  . 

The  key  to  the  solution  of  the  problem  lies  in  Eq.  (4c).  By 

adjusting  the  end  point  energy  of  the  irradiation  to  successively 

lower  values,  the  F  (E.)  can  be  made  zero  for  all  but  the  lowest 

m  l 

energy  gateway,  E1  .  Then,  there  is  only  one  term  on  the  right. 
Further  adjustment  of  the  end  point  will  identify  that  gateway 
energy  by  turning  off  the  yield,  as  seen  in  the  example  shown  in 
Fig.  2.  There,  E^  is  clearly  indicated  by  the  intercept. 
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U,  r 


END  POINT  ENERGY  (MeV) 


■g-jre  2:  Nuclear  fluorescence  from  ^ 1  ^mln  normalized  to 
*XBr  pumped  by  single  pulses  of  bremsstrahlung  having  end 
point  energies  shown.  The  linear  fit  to  the  data  intercepts 
the  axis  at  a  gateway  energy  of  1.078  MeV. 


Once  the  lowest  gateway  energy  is  determined,  Eq.  (4c)  can 
solved  for  the  first  cross  section, 

Q,  =  S  /  F  (  E.  )  ,  (5) 

i  m  ml 

ere  the  free  index,  m  implies  that  Q-  can  be  determined  as 


an  average  from  a  number  of  measurements.  It  is  not  rigorously 
true  that  spectral  content  can  never  be  measured  for  any 
bremss trahlung .  If  intense  pulsed  sources  can  be  fired  with 
reasonable  cost  and  repetition  rate,  they  can  be  observed  with  a 
geometry  insuring  that  only  one  photon  reaches  a  detector  per 
pulse.  In  that  case  the  energy  of  that  photon  can  be  determined 
and  logged  before  the  next  pulse  and  in  this  way  the  spectrum  of 
the  bremsstrahlung  can  be  collected  as  an  experimental 
probability  of  detecting  a  photon  of  a  given  energy.  From  the 
geometry  the  total  number  of  photons  can  be  determined  at  a 
particular  energy  for  a  given  duration  of  the  m-th  exposure 
producing  the  measured  number  of  activations,  S  .  Then,  of 

course,  Eq.  (5)  can  be  used  to  extract  a  value  of  which  can 

subsequently  be  used  to  determine  fluences  at  E.^  from  any  other 
sources  by  measuring  only  the  S  they  produce. 

Once  the  first  Q1  is  found,  Eq.  (4b)  can  be  used  to 

determine  the  next  with  measurements  made  with  a  higher  end 
point  energy, 

°2  ■  <sm  -  °lW>  '  W  ‘  <61 

In  actual  practice  a  great  amount  of  care  is  necessary  to  insure 

that  changes  of  end  point  energy  do  not  open  more  than  one 

gateway  at  a  time  and  to  collect  sufficient  statistics  for  the 

F  (E.)  so  that  the  differences  implied  by  Eq.  (6)  do  not  degrade 
m  l 

the  determinations  of  the  Q^. 

It  appears  that  only  variable  energy  linacs  have  the 
combinations  of  parameters  needed  to  meet  these  requirements  and 
one  must  be  procured  in  subsequent  stages  of  this  development. 
Once  the  equivalent  of  Table  I  is  produced  in  this  way  for  the 
range  of  energies  from  1.4  to  6.0  MeV,  the  continuum  from  any 
source  can  be  measured  directly  with  a  simple  set  of  foils  and  a 


counting  system. 


Ill  Phase  I  Accomplishments 


Our  Phase  I  modeling  has  shown  that  gateway  energies  for  the 
activation  of  1 9  stable  isotopes  could  be  measured  with  a 
variable  energy  linac  which  will  be  needed  for  the  next  phase  of 
this  research.  The  scaling  of  preliminary  measurements  with  a 
fixed  energy  linac  has  proven  that  sensitivities  of  all  19 
materials  are  sufficient  for  use  in  sampling  the  spectral  content 
of  a  single  shot  from  a  nuclear  simulator  of  the  scale  of 
DNA/ AURORA,  as  well  as  the  spectrum  from  a  linac  of  the  type  used 
for  medical  radiation  therapy  during  the  few  minutes  period  each 
day  set  aside  for  calibration.  The  viability  of  a  commercial 
instrument  or  service  seems  very  clear,  once  the  nuclear 
properties  of  the  19  are  standardized  in  subsequent  work. 

The  principal  risk  is  that  all  isotopes  might  conceivably 
sample  the  same  slice  of  the  spectrum  and  so  be  unable  to  record 
spectral  intensities  at  energies  other  than  at  that  single  value. 
Were  there  only  two  or  three  usable  isotopes,  such  a  possibility 
might  be  cause  for  concern.  However,  the  chances  that  19 
isotopes  would  have  indistinguishable  absorption  spectra  is  too 
remote  for  serious  consideration.  Such  concern  is  more  properly 
reduced  to  the  status  of  a  problem;  a  problem  in  optimizing  the 
best  spectral  coverage  with  the  most  economical  target  package. 
That  is  part  of  the  work  needed  in  a  future  effort. 

The  details  of  the  accomplishments  of  Phase  I  upon  which 
rest  the  strongly  positive  assessment  of  feasibility  are  reviewed 
below. 


1  3 


III.1  The  19  sampling  isotopes. 

Appearing  in  Table  II  is  a  proprietary  listing  of  the  19 
stable  isotopes  known  to  activate  when  irradiated  by  the  fluence 
accumulated  from  a  6  MeV  Varian  CLINAC  1800  in  a  few  minutes. 
Shown  there  is  the  isomer  and  its  abundance  in  naturally 
occurring  material  in  the  first  two  columns.  In  the  next  columns 
are  the  energies  of  the  isomeric  levels  and  their  lifetimes  for 
storing  the  samples  of  the  irradiating  spectrum  to  which  they  are 
sensitive.  Then  follow  statistical  quantities  for  use  in 
Eq.  (3a)  with  the  identification.  Ig  =  and  Ig  =  Jg.  Finally, 
in  the  last  column  of  Table  II  is  the  most  important  quantity, 
the  integrated  cross  section,  q  of  Eq.  (3b). 

Since  the  energy,  E^  of  the  gateway  for  the  activation  of  a 
particular  isomer  is  not  known,  the  individual  cross  sections 
cannot  be  expressed  as  the  q^  for  a  given  E^ .  That  is  the  problem 
onto  which  the  next  phase  of  work  must  be  focused.  Rather,  in 
Table  II  the  artifice  has  been  employed  of  computing  the  value  of 
q  which  would  have  been  required  to  produce  the  observed 
activation  through  a  gateway  of  2.125  MeV  which  is  most  probably 
the  lowest  value  at  which  such  a  large  activation  could  occur. 

The  cross  sections  of  Table  II  were  determined  from  the 
absolute  activations  resulting  from  irradiations  with  the 
spectrum  shown  in  Table  III.  As  can  be  seen,  it  has  been 
normalized  to  unit  spectral  intensity  at  2.125  MeV.  Had  the 
gateway,  E^  in  a  particular  case  opened  at  a  higher  energy,  then 
the  observed  amount  of  product  would  have  resulted  from  a 
spectral  intensity  lower  than  assumed  and  the  actual  value  of 
cross  section  would  be  larger  than  what  is  shown  in  Table  II  by 


Table  II 


Summary  of  the  properties  of  isomers  known  to  activate 
through  giant  gateways  opening  in  the  range  of  energies  between 
2  MeV  and  6  MeV.  Values  for  q,  the  integrated  cross  section  of 
Eq.  (3b)  are  taken  from  Ref. 16.  Since  the  gateways  are  unknown, 
the  cross  sections  have  been  reported  as  values  pertaining  to  a 
standard,  arbitrary  reference  energy  of  2.125  MeV. 
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Table  III 


Ratios  of  the  photon  intensities  at  energies  E  to  that  at 
2.125  MeV  for  the  spectrum  used  in  the  measurements  of  Table  II. 


Energy  r  =  4>(  E)  /4>(  2  .  1  25  ) 

(MeV) 


0.125 

0.015 

0.375 

0.500 

0.625 

1  .731 

0.875 

1  .769 

1.125 

1  .  654 

1  .  375 

1  .538 

1  .  625 

1  .538 

1  .875 

0.962 

2.125 

1  .  000 

2.375 

0.731 

2.625 

0.538 

2.875 

0.500 

3.125 

0. 61  5 

3.375 

0.354 

3.625 

0.500 

3.875 

0.385 

4.125 

0.346 

4.375 

0.204 

4.625 

0.162 

4.875 

0.188 

5.125 

0.065 

5.375 

0.1  08 

5.625 

0.050 

5.875 

0.065 
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the  compensating  amount.  For  example,  if  a  particular  were 
found  to  be  2.875  MeV  rather  than  the  arbitrary  reference  value 
of  2.125  then  Table  III  shows  that  the  actual  integrated  cross 
section,  q^  should  be  twice  the  value  of  q  shown  in  Table  II. 
Generally,  if  there  is  single  gateway  at  ,  then, 

q±  =  q/r(E.)  .  (7) 

The  situation  with  multiple  gateways  is  more  complex  and  the 
equivalent  expression  is  found  from  the  relation, 

q  =  r(E1)q1  +  r(E2)q2  +  ....  (8) 

The  decomposition  of  q  into  the  individual  qi  can  only  be  done 
experimentally,  as  expected  in  a  future  phase  of  effort.  However, 
nuclear  systematics  provide  guidance  sufficient  to  show  the  clear 
feasibility  of  such  experiments.  Moreover,  they  serve  to 
identify  seven  candidates  for  the  role  of  primary  norm. 

III. 2  Nuclear  Systematics  and  the  Primary  Norm 

Completed  during  this  Phase  I  of  research  has  been  a  study 
of  the  systematics  of  the  nuclear  properties  contributing  to  the 
occurrence  of  giant  gateways  for  activation  such  as  those  shown 
in  Table  II.  Figure  3  shows  a  plot  of  the  integrated  cross 
section,  q  as  a  function  of  the  mass  number,  A  of  the  nucleus 
being  activated. 

Of  great  importance  in  Fig.  3  are  the  two  triangular  data 

points  which  record  the  only  two  cases  for  which  the  E^  have  been 

determined  for  excitation  of  isomers  through  giant  gateways. 

17  1  8  87  1  97 

Those  two  isotopes  '  were  Sr  and  Au.  For  consistency, 
the  integrated  cross  sections  in  both  cases  have  been  "corrected" 
back  to  2.125  MeV  for  plotting  through  use  of  the  inverse  of 


19*/ 

Eq.  (7).  The  data  point  for  'Au  is  shown  as  a  part  of  the 

1  8 

total  q  of  Eq.  (8)  that  was  separated  by  an  ad  hoc  procedure 
applicable  only  to  that  material. 


MASS  NUMBER  (A) 


Figure  3:  Plot  of  the  integrated  cross  sections,  q  as 

functions  of  the  mass  number  of  the  nuclei  listed  in  Table 

II.  Triangles  plot  data  points  for  individual  q.  where 

8^7 

values  are  available  from  Refs.  17  and  18  for  Sr  and 
197Au,  respectively. 

Even  a  cursory  inspection  shows  several  trends  worthy  of 
note.  With  the  exception  of  the  anomalously  large  value  for 
77Se,  the  other  18  isotopes  seem  to  group  into  three  islands,  the 
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peaks  of  which  appear  to  scale  up  by  about  a  factor  of  five  from 
the  one  immediately  preceding.  Countering  this  general  trend  is 
a  distribution  of  much  smaller  values  which  seem  to  vary  little 
over  the  whole  span  of  masses. 

Since  the  funneling  process  of  Fig.  1  should  not  work  so 
well  for  purely  single  particle  states  of  nuclear  excitation,  it 
has  been  hypothesized  that  collective  excitation  of  several  core 
nucleons  could  break  symmetries,  tending  to  hinder  the 
transitions  of  Fig.  1  because  of  selection  rules  on  projections 
of  angular  momenta.  At  the  same  time,  the  collective  excitation 
of  several  particles  would  tend  to  enlarge  the  magnitudes  of  the 
matrix  elements  coupling  the  nucleus  to  the  radiation  field  so 
that  transitions  would  occur  with  the  same  multipolarities  as 
indicated  by  the  AJ  for  the  transition. 

The  general  scaling  of  electromagnetic  transition  rates  when 

single  particles  mediate  the  transitions  is  given  by  a  number  of 
21 

authors  in  forms  equivalent  to, 

W ( QUAD )  =  CqA4/3E5  ,  (9) 

for  quadrupole  transitions,  and 

W(OCT)  =  CqA2E7  ,  (10) 

for  octupole  transitions.  While  it  is  extremely  speculative  to 
attempt  to  fit  cause  to  the  experimental  effects  recorded  in 
Table  II,  there  is  no  theoretical  guidance,  so  perhaps  some 
empiricism  is  justified.  There  is  nothing  else  possible  until 
further  work  is  completed. 

20 

It  has  been  reported  that  collective  octupole  transitions 
should  be  much  stronger  than  quadrupoles  because  the  symmetries 
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of  such  vibrations  tend  to  release  the  otherwise  hindered 

electric  dipole  components  of  a  transition.  In  that  case  we 

might  reasonably  identify  the  smaller  distribution  of  cross 

sections  with  the  less  effective  collective  quadrupole 

vibrations.  Then,  using  tne  known  values  shown  in  Fig.  (3)  for 

87  197 

tne  two  oat  a  points,  r(E^)q^  for  Sr  and  Au  allows  for  the 
determination  of  C^E3  in  Eq .  (9).  The  resulting  fit  is  shown  in 
Fig.  4  where  the  dotted  curve  plots  the  value, 

v:  =  1  000(A/87)4/3  ,  (11  ) 

-29  2 

in  the  usual  units  of  1 0  cm  keV. 

A  magnification  of  the  relevant  portions  of  Fig.  4  is  show’n 
in  Fig.  5.  The  agreement  to  the  expression  of  Eq.  (11)  of  the 
data  for  the  seven  isomers  specifically  identified  is  much  better 
than  could  be  reasonably  expected,  were  there  not  some  merit  in 
the  approximations  being  employed. 

As  a  working  hypothesis  for  further  analysis,  it  was  assumed 

that  the  seven  elements  summarized  in  Table  IV  were  activated 

through  a  single  gateway  lying  in  the  quadrupole  region  between 

about  1.5  and  3.0  MeV.  The  much  greater  total  cross  sections  of 

the  other  isotopes  were  then  attributed  to  the  contribution  of  a 

larger  gateway  opening  at  higher  energies.  This  was  known  to  be 

197 

the  case  for  the  only  system  studied  in  detail  previously,  Au. 
For  use  in  our  computations,  the  q^  for  the  seven  isotopes  most 
likely  to  be  single  gateway  materials  were  obtained  from 

qi  =  1 000(A/87)4/3/r (Ei )  .  (12) 

The  role  of  such  single  gateway  isotopes  as  normalizing  standards 
is  extremely  important  in  the  measurement  of  spectral 
intensities,  and  it  is  encouraging  to  reflect  that  success  in 
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implementing  the  XAN  process  over  the  lower  energy  range,  0.2-1. 5 

79 

MeV  required  only  one  such  selective  nuclei,  Br  over  that 
range.  Here,  for  this  case  there  are  seven  potential  standards. 
Any  one  of  them  would  be  viable  as  a  primary  norm.  Probabilities 
are  very  high  that  one  can  be  proven  experimentally  while  the 
others  are  used  to  provide  convenient  spectral  detail. 
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Figure  4:  Plot  of  the  integrated  cross  sections,  q  as 
functions  of  the  mass  number  of  the  nuclei  listed  in  Table 
II.  The  dotted  curve  shows  the  scaling  of  transitions 
mediated  by  a  giant  quadrupole  oscillation. 
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Figure  5:  Magnification  of  a  portion  of  Fig.  4  plotting 
values  of  the  integrated  cross  sections,  q  from  Table  II  as 
functions  of  mass  number.  Isotopes  are  identified  whose 
cross  sections  seem  to  agree  with  the  approximation  of 
Eq.  (11)  shown  by  the  dotted  curve. 


Table  IV 


Isotopes  identified  in  Phase  I  studies  as  the  roost  likely  to 
have  single  gateways  for  activation  in  the  2-6  MeV  range.  These 
are  the  trial  standards  for  subsequent  validation.  Recorded  in 
the  rightmost  column  are  the  integrated  cross  sections,  q 
obtained  from  the  smoothing  approximation  for  quadrupole 
transitions,  Eq.  (12)  normalized  to  2.125  MeV  for  convenience  by 
Eq.  (7). 
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II. 3  Measurement  Strategies 


Even  an  a  priori  choice  of  a  primary  norm  from  the  group  of 
seven  candidates  would  be  dominated  by  what  is  now  a  free  choice 
cf  the  functional  cycle.  There  are  two  classes  of  operation  to 
consider.  The  most  versatile  is  a  "pop-in"  strategy  which  is 
also  the  most  expensive  to  engineer  as  a  subsequent  product  or 
service.  There  the  transfer  from  the  point  of  irradiation  to  the 
place  of  measurement  is  done  pneumatically  or  mechanically.  With 
such  a  system,  the  time  lost  in  transfer  can  be  as  little  as  a 
second  so  that  even  the  isotopes  with  the  shortest  storage  times 
car.  retain  enough  activation  for  accurate  read-out.  However,  the 
complexity  cf  such  a  system,  is  so  much  greater  than  for  the 


alternative,  any 

strategy 

for 

commercialization 

would 

require 

first  preference 

be 

giver 

to 

the  other 

technique.  Because 

of 

this,  the  pnnci 

pal 

emph 

asis 

in  Phase 

I  was 

placed 

on 

the 

alternative  cycle  with  the  intention  to  switch  to  the  "pop-in" 
approach  if  it  became  clear  that  either  accuracy  or  resolution 
was  being  unduly  compromised.  That  necessity  has  not  yet  been 
perceived . 

The  preferred  strategy  for  the  implementation  of  XAN  is  a 
passive  one  in  which  a  target  package  is  placed  in  the  area  to  be 
irradiated  by  an  operator  or  technician.  Allowing  operator 
choice  permits  the  precise  point  of  medical  treatment  or 
simulated  blast  impact  to  be  instrumented  with  the  ultimate 
flexibility.  After  irradiation  the  operator  returns  to  detach 
the  target  and  insert  it  into  a  counting  chamber  which  is 
properly  shielded  to  insure  a  low  background.  An  answer  would  be 
desirable  within  a  time  comparable  to  that  for  irradiation, 
transport  and  handling,  and  that  implies  only  a  part  of  the 
activation  will  be  counted  during  that  arbitrary  standard 
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interval.  Each  step  of  treatment  is  associated  with  a  loss  of 
activation ,  and  hence,  a  reduced  sensitivity  that  will  change  the 
figures  of  merit  associated  with  each  of  the  candidate  norms  of 
Table  IV  from  being  simply  dependent  upon  the  value  of  q(2.125). 

For  the  m.odeling  of  the  efficiencies  for  the  use  of  any  of 
the  candidate  norms,  it  is  necessary  to  obtain  an  effective  value 
of  integrated  cross  section,  c'  that  has  been  corrected  for  the 
varying  amounts  of  loss  of  activation  imposed  by  the  particular 
operating  cycle  planned, 


c'eff  *  abcs  ,  (13) 

where  b  is  the  loss  of  activation  during  transfer, 

b  =  exp(-wtt)  ,  (14a) 

where  c  is  the  fraction  of  the  activation  which  fluoresces  during 
the  counting  period, 

c  =  l-exp(-wt  )  ,  (14b) 

where  w  is  the  storage  lifetime  of  the  activation  and  and 
t  are  the  times  for  transfer  and  counting,  respectively.  The 
correction  factor,  a  for  the  decay  of  the  activation  during  the 
finite  duration,  t  of  the  irradiation  is  unity  for  the  case  of 
impulsive  excitation,  but  otherwise  is 

a  =  (wtr )-1 ( 1 -exp(-wtr ) )  ,  (14b) 

where  a  has  been  written  to  correct  activation  per  unit  dose  (as 
opposed  to  dose  rate). 


Choosing  the  modeling  parameters  for  validating  medical 
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machines, 


irradiation  time,  t  =10  min. 

r 

transfer  time,  t ^  =  2  min. 

counting  time,  t  =10  mm., 


(15a) 


1  5c) 


values  which  create  little  operator  stress,  the  resulting  figures 
of  merit  for  the  candidate  norms  can  be  computed  from  Eq.  (13). 
However,  another  important  aspect  warrants  inclusion  in  this  type 
cf  effective  cross  section  and  that  is  the  fractional  abundance, 
f  of  the  isotope  so  that 

a  =  fq’  r,  =  abcfq  .  (16) 

_el l  eft 

This  quantity  more  nearly  approximates  the  useful  activation  per 
input  ccse  per  unit  amount  of  target  material.  It  is  tabulated 
m  Tarle  V  for  the  seven  candidate  norms.  The  comparative  values 
for  impulsive  excitation,  for  which  a  =  1,  are  also  shown. 

III. 4  Preferred  Norrr.s 

It  is  clear  from  Table  V  that  two  of  the  seven  isotopes  are 

preferred  choices  for  use  as  primary  norms,  Ba  and  Hg. 

Even  with  the  rather  relaxed,  passive  strategy  for  measurement 

the  effective  cross  sections  for  activation  remain  in  the  10-100 

(x10  cm  keV)  region  we  have  characterized  as  describing  giant 

quadrupole  transitions.  Both  considerably  exceed  the  value  of 

7  9 
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6.2  which  was  sufficient  for 


Br  to  serve  as  a  norm  in  the 


energy  range  below  1.5  MeV.  However,  they  differ  from  each  other 
in  the  details  of  usage. 


The  mercury  isotope,  Hg  appears  to  be  the  ideal  norm.  No 
other  isotope  of  mercury  activates  and  the  time  for  decay  of  the 
signal  is  long  enough  that  several  successive  read-out  cycles 


would  be  possible  in  the  event  operator  error  were  suspected. 
Conversely,  it  is  not  so  long  that  the  sample  could  not  be  reused 
the  next  day  without  the  need  for  making  a  baseline  measurement 
of  residual  activation  before  sampling.  If  more  frequent 
validation  of  the  source  were  desired,  mercury  is  inexpensive 
enough  that  several  samples  could  be  stocked  to  avoid  the 
nuisance  of  counting  residual  activation. 

Table  V 

The  effective  values  of  integrated  cross  sections  for  the 
isotopes  possibly  useful  as  primary  norms  which  relate  signal 
read-out  to  applied  dose  in  a  passive  sampling  geometry  creating 
low  operator  stress. 


1 s  ctcpe 

Environment 

q  Medical  Therapy 

el-L  _29  2 

(X"iC  'err.  keV) 

q  Simulation 

Gir  -29  2 

(XI  0  cm,  keV) 

0.00 

0.00 

e7Sr 

2.76 

2.81 

65v 

0.22 

5.80 

,17sr. 

0.04 

0.04 

137Ba 

39.90 

116.5 

1  6  3t, 

0.00 

0.00 

1  9  9U 

Hg 

68.00 

73.67 
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The  1  J  '  Ba  is  even  more  favorable  in  these  regards,  and  its 
desirability  as  a  primary  norm  is  lessened  only  by  the  activation 

I-35 

of  the  ”  Ba  component  of  natural  barium  present  in  comparable 

abundance.  This  latter  constituent  has  a  half  life  over  a  day 

and  will  gradually  accumulate  activation  which  will  contribute 

some  noise  to  the  counting  environment.  However,  since  the  more 

comprehensive  measurement  strategies  anticipate  the  use  of  1 ^Ba 

as  a  multiple  gateway  material,  the  recycling  of  this  isotope 

appears  inevitable.  As  long  as  the  bookkeeping  is  done  to  record 

residual  activation,  then  "*  "^Ba  becomes  almost  as  attractive  a 

1  99 

primary  norm,  as  Hg. 

The  conclusion  of  this  part  of  the  Phase  I  analysis  is: 

1  99 

Primary  Norm,  -  Hg 

1  ,7 

Secondary  Norm  -  Ba . 

III. 5  Quadrupole  and  Octupole  Transitions 

For  the  12  isotopes  appearing  in  Table  II  that  do  not  appear 

in  Table  IV,  the  indication  is  that  their  much  larger  value  of 

integrated  cross  section  is  contributed  by  the  opening  of  a  giant 

octupole  resonance  at  higher  energy.  There  may  or  may  not  also 

be  a  quadrupole  resonance,  so  that  these  are  most  probably 

1  97 

multiple  gateway  materials.  That  is  the  case  of  Au,  the  only 

isotope  for  which  such  gateways  have  been  approximately 

,  «.  .  1  8 
located . 

The  relative  size  of  octupole  cross  sections  to  the 
quadrupoles  is  so  great  that  plots  of  total  q  can  be  compared  to 
parameterizations  such  as  Eq.  (10)  for  octupole  transitions  with 
little  error,  even  if  the  quadrupole  contributions  are  not 
removed.  In  Fig.  6  the  approximations  of  Eq.  (10)  are  plotted  by 
solid  curves  for  two  cases  that  might  reasonably  bound  the 
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Figure  6:  Plot  of  the  integrated  cross  sections,  q  as 
functions  of  the  mass  numbers  of  the  nuclei  listed  in  Table 
II.  The  solid  curve  shows  the  scaling  of  transitions 
mediated  by  giant  octupole  oscillations.  Energies  of  the 
octupole  gateways  are  given  in  terms  of  the  energies  of 
quadrupole  transitions  for  the  same  isotope. 


majority  of  isotopes.  Recalling  that  data  in  Fig.  6  are  scaled 
down  to  2.125  MeV  by  Eq .  (8),  then  only  the  ratio  of  the  EQ  for 
the  octupole  resonance  to  the  E^  for  the  quadrupole  resonance 
need  be  raised  to  the  seventh  power  in  Eq.  (10)  so  that 

=  500(A/87)2(  E_/E_. )  7  , 


(17) 


where  the  value  of  has  been  chosen  to  be  500  to  fit  the  data 
near  A  =  200.  As  seen  in  the  figure,  the  octupole  gateways  would 
seer,  to  lie  between  1.25  E^  and  1.5  E^.  Again,  the  level  of 
agreement  is  surprising  just  as  in  the  case  of  fitting  the  single 
quadrupoie  gateways  in  Fig.  5.  However,  in  this  case,  the 
exponent  of  the  dependence  on  energy  is  so  great  that  actual 
values  are  to  be  preferred  over  those  given  by  Eq .  (17). 

III. 6  Measurement  Modeling 

An  important  concern  is  whether  there  are  enough  of  the 
isotopes  with  octupole  transitions  and  storage  lifetimes  long 
enough  to  be  compatible  with  the  passive  measurement  strategy 
mentioned  above.  Table  VI  summarizes  the  effective  cross 
sections  from.  Eq.  (16)  for  those  isotopes  of  Table  II  having 
octupole  transitions  not  reduced  to  zero  by  the  losses  in 
transfer,  and  by  the  brevity  of  the  counting  interval. 

Recalling  that  Table  I  shows  cross  sections  as  small  as  0.2 
to  have  been  useful  in  measuring  spectral  intensities  at  energies 
below  1.5  MeV,  it  can  be  seen  from  Table  VI  that  there  are  seven 
usable  octupole  transitions,  even  for  the  most  convenient 
strategy  of  measurement.  The  full  cross  sections  of  10,000’s  of 
units  would  be  available  in  a  "pop-in"  system,  but  as  yet, 
analysis  does  not  mandate  the  extra  difficulties.  It  seems 
preferable  to  exploit  the  large  size  of  the  octupole  transitions 
to  create  enough  activation  that  some  can  be  wasted  by  the 
simpler  procedures  for  handling  the  samples  that  would  be  most 
useful  in  a  practical  device. 

Since  each  of  the  materials  of  Table  VI  is  presumed  to  have 
at  lease  two  energies  Ei  at  which  spectral  intensities  are 
sampled,  there  are  14  points  from  the  isotopes  of  that  list  and 
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another  four  from  Table  V.  In  all,  there  are  expected  to  be  18 
points  at  which  the  intensities  are  recorded. 

Table  VI 

The  effective  values  of  integrated  cross  sections  for  the 
isotopes  apparently  having  octupole  transitions.  These  quantities 
relate  signal  read-out  to  applied  dose  in  a  passive  sampling 
geometry  creating  low  operator  stress. 


Isotope 


Environment 


c  ,,,  Medical  Therapy 
'err  _pg  2 

(XI 0  cm  keV) 


qeff ' 
(XI  0 


Simulation 

-29  2. 

cm  keV) 


77Se 

0.21 

5.00 

n1Cd 

50.9 

54.6 

113- 

_n 

14.8 

14.8 

115ln 

1  64 

1  64 

135Ba 

1  .7 

1.7 

1  7  '  H  f 

2.2 

49.2 

195Pt 

8.5 

8.5 

If  we  imagine  the  sampling  points  arranged  in  an  ascending 
series  of  activation  energies,  they  will  comprise  a  set  of 
energies,  where  i  runs  from  1  to  18.  They  will  be  distributed 
over  a  range  from  2  MeV  to  some  maximum  assumed  to  be  6  MeV.  The 
differences 


where  EQ  =  2  MeV,  then  comprise  a  set  of  random  variables  whose 
sum  is  conserved  at  4  MeV.  The  average  is  X  =  0.222  MeV.  These 
intervals  X.^  and  sampling  points  follow  the  same  statistics  as 
particles  in  an  ideal  gas  and  their  kinetic  energies.  It  can  be 
shown  that  the  relative  probability  of  a  particular  X.^  is 

Pi  =  exp  (-Xi/X)  ,  (19) 

and  the  normalized  probability,  should  be, 

dPi  =  exp(-Xi/X)d(Xi/X)  .  (20) 

The  integration  of  Eq.  (20),  then  shows  there  to  be  a  proability 
of  only  11%  that  there  will  be  a  gap  larger  than  0.5  MeV  that  is 
not  sampled  by  one  of  the  18  transitions  of  the  isotopes  listed 
in  Tables  V  and  VI. 

This  seems  a  reasonable  likelihood  of  an  adequate  resolution 
for  the  measuring  process  being  modeled. 

The  actual  computational  procedure  to  be  followed  for 
unfolding  the  spectrum  from  the  set  of  activations  it  produces  in 
the  11  isotopes  used  for  sampling  will  depend  upon  the  particular 
energies  sampled  by  the  single  gateway  materials.  Depending  upon 
their  arrangement  in  the  interval  of  energies  from  2  to  6  MeV, 
the  unfolding  process  will  be  either  interpolative  or 
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extrapolative,  or  more  likely,  a  combination  of  these  techniques. 
The  following  example  is  illustrative. 


The  spectrum  of  no  source  is  well  enough  known  to  be  used 
without  further  arbitrary  assumptions.  However,  at  least  the 
distribution  of  energies  sampled  is  well  known  from  Table  I  over 
the  range  from  0.2  to  1.5  MeV,  so  this  range  is  the  one 
preferable  for  an  example.  For  these  energies  the  linac  spectrum 
of  Table  II  can  be  simply  scaled  along  the  abscissa  to  fit  an 
endpoint  of  1.5  MeV,  as  shown  in  Fig.  7. 


PHOTON  ENERGY  (  MeV  ) 


Figure  7:  Plot  of  a  hypothetical  linac  spectrum  obtained  by 
scaling  the  data  of  Table  II  to  the  energy  interval  shown. 
The  solid  curve  shows  the  arbitrary  fit  of  a  smooth  curve  to 
this  data. 


33 


To  avoid  introducing  a  favorable  bias  which  might  occur 
through  use  of  a  power  series,  an  arbitrary  curve  was  fit  by  hand 
to  the  data  of  Fig.  7  as  shown  by  the  solid  curve.  Had  a  power 
series  been  used,  then  the  same  order  of  series  would  have 
allowed  the  subsequent  unfolding  to  be  exact,  and  so  minimized 
the  intrinsic  difficulties  of  the  problem.  In  contrast,  a  fit  by 
hand  gives  an  unknown  mathematical  form. 

Neglecting  the  minor  contribution  of  the  gateway  at  250  keV 
77 

in  Se,  the  matrix  for  the  computation  of  Eq.  (1)  can  be  written 

as  shown  in  Table  VII.  The  activation  of  a  particular  isotope, 
77 

Se  for  example,  can  be  seen  to  be 

51  =  Q1F(E1 )+Q3F(E3)+Q4F(E4)  (21a) 

while  the  other  activations  occur  through  single  gateways,  namely 

52  =  Q2F(E2)  ,  (21b) 

53  =  Q5F(E5)  ,  (21c) 

54  =  Q6F(E6)  .  ( 21 d ) 

This  set  of  equations  is  undrrdetermined  and  cannot  be  simply 
inverted  because  the  spectrum,  F(E)  is  being  sampled  at  six  Ei  to 
yield  only  four  measurements. 

The  importance  of  the  single  gateway  materials  is  readily 
apparent  from  Eqs.  ( 21 b ) - ( 21 d ) .  •  These  can  be  immediately 
inverted  to  give  F(E2),  F(Eg)  and  F( Eg ) .  The  basic  assumption 
necessary  to  obtain  an  approximate  solution  for  F(E^),  F(E3)  and 
F(E4)  is  that  the  function  F  be  continuous  with  continuous 
derivatives.  In  that  case  F(E3)  and  F(E4)  can  be  obtained  by 
interpolating  between  the  points  obtained  from  the  single  gateway 
isotopes.  Using  a  polynomial  fit  to  F(E2),  F(Eg)  and  F(Eg)  gives 
reasonable  estimates  for  F(E-j)  and  F(E4)  as  shown  in  Fig.  8, 
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where  the  results  of  this  unfolding  process  are  shown  together 
with  the  original  data  for  the  linac  spectrum. 

The  measurement  at  FfE^)  is  pathologically  difficult  and 
makes  a  good  example.  The  extrapolation  of  polynomial  fits  is  a 
notoriously  poor  procedure.  Since  extrapolation  is  needed  for 
only  one  point,  an  alternative  is  to  solve  Eq.  (21a)  after 
substituting  the  values  interpolated  for  F(E^)  and  F(E^). 
Unfortunately,  so  this  process  requires  the  result  be 
obtained  as  a  small  residue  between  large  measurements. 
Nevertheless,  even  in  this  pathological  case,  the  final  value 
obtained  for  F(E.j)  is  not  too  poor,  as  can  be  seen  in  Fig.  B. 
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Figure  8:  Plot  of  the  hypothetical  linac  spectrum  shown  in 
Fig.  7,  together  with  the  values  of  the  six  measurements 
unfolded  from  simulated  activations. 
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This  example  considered  in  Phase  I  is  a  worst  case,  and  it 
is  expected  that  the  distribution  of  gateways  in  the  1. 5-6.0  MeV 
range  will  be  much  more  favorable.  As  discussed  in  previous 
sections,  it  is  the  selection  of  isotopes  of  Table  IV  that  were 
concluded  to  have  single  gateways.  Contributed  by  the  quadrupole 
transitions  at  lower  energies,  the  activations  should  be  much 
smaller  than  what  results  from  the  octupole  transitions  at  higher 
energies.  Even  if  there  would  be  only  multiple  gateways  in  the 
octupole  region,  extrapolation  in  this  case  would  consist  of  the 
removal  of  small  fitted  contributions  from  large  total 
activations  instead  of  the  inverse  situation  encountered  in  the 
example  reported  above.  In  contrast  to  that  pathological  case, 
the  situation  for  the  energy  range  2.0  -  6.0  MeV  has  been  shown 
by  the  Phase  I  work  to  be  far  more  tractable.  Nevertheless,  even 
at  the  more  difficult  lower  energies,  the  use  of  XAN  to  measure 
spectral  intensities  has  already  been  proven  in  numerous 
examples . 

The  modeling  of  specific  hypothetical  cases  was  facilitated 
by  the  development  of  two  software  packages: 

1 )  ENFOLD  to  generate  synthetic  activation  data,  and 

2)  UNFOLD  to  invert  it  in  order  to  obtain  the  irradiating 

spectrum  originally  assumed. 

Run  against  each  other  they  were  able  to  demonstrate  that  for 
likely  values  for  the  unknown  nuclear  parameters  described  above 
unfolded  spectra  could  be  realized  with  acceptable  statistics 
when  used  in  the  targeted  measurement  strategy.  Both  programs 
were  completed  during  this  phase  of  research  and  documentation  is 
reproduced  in  Appendix  A.  As  can  be  seen  there,  the  input 
spectrum  F  (E.)  can  be  entered  either  as  a  table  of  values  for 

ra  l 

pathological  cases  or  may  be  selected  with  a  simple  keypress  to 
correspond  to  a  Varian  CLINAC  (6MeV),  to  AURORA,  or  to  our 
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estimates  for  the  output  of  a  variable  energy  linac  with  end 
points  being  entered  in  an  interactive  mode. 


Careful  attention  was  given  to  the  implementation  of  the 
correct  statistical  noise  distribution.  Since  small  activations 
are  economical  in  terms  of  the  amount  of  material  required  and 
the  duration  (or  proximity)  of  the  calibrating  period,  it  was 
important  to  use  Poisson  distributions  of  counting  rates.  Once  an 
activation  was  calculated  in  a  particular  circumstance,  a  random 
amount  of  noise  conforming  to  Poisson  statistics  was  added  by  the 
ENFOLD  program  to  model  the  real  uncertainty  in  the  results  of 
counting  the  decay  of  a  single  channel  of  activation. 

In  Appendix  A  can  be  seen  the  details  of  the  use  of  these 
programs,  together  with  results  in  illustrative  cases.  The 
unfolding  of  the  hypothetical  spectrum  can  be  seen  to  be 
remarkably  successful  in  the  various  practical  cases  of  interest. 


IV.  Summary  of  Accomplishments 


As  detailed  above,  the  achievements  realized  in  Phase  I  of 
this  research  prove  the  clear  feasibility  of  subsequent  work.  A 
primary  norm  having  a  single  gateway  for  activation  has  been 
found  for  the  2.0  to  6.0  MeV  range  of  Photon  energies,  and  eleven 
isotopes  have  been  identified  as  able  to  sample  the  spectrum  at 
18  points.  The  modeling  of  a  preferred  measurement  strategy  has 
shown  that  atarget  composed  of  these  materials  could  be  manually 
positioned  by  an  operator,  exposed,  and  then  transferred  by  the 
operator  to  a  measurement  station  located  in  reasonable 
proximity.  Counting  times  of  only  10  minutes  and  exposures  of 
either  10  minutes  for  a  medical  linac  or  one  shot  from  a  nuclear 
simulator  would  give  adequate  activation  for  analysis.  The 


modeling  strategies  further  showed  the  probabilities  were  only 
10%  that  a  gap  as  large  as  0.5  MeV  would  be  left  unmeasured  by 
the  selection  of  isotopes  prescribed,  and  that  there  was  far  less 
opportunity  for  a  pathological  accident  in  the  distribution  of 
gateways  than  in  the  case  of  measurement  below  1.5  MeV,  already 
reduced  to  practice. 

The  only  remaining  obstacle  to  the  actual  realization  of  an 
XAN  procedure  for  validating  bremsstrahlung  spectra  in  the  2.0  to 
6.0  MeV  region  lies  in  the  acquisition,  installation  and 
operation  of  a  variable  energy  linac  with  which  the  gateways  for 
sampling  spectral  intensities  can  be  pinpointed. 
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Appendix  1 


Review  of  the  ENFOLD  and  UNFOLD  programs  for  generating  and 
deconvolving  synthetic  data. 

A  maior  project  goal  for  Phase  I  was  met  by  the  development 
of  two  programs  which  allowed  the  simulation  of  a  spectral 
calibration.  The  first,  called  ENFOLD,  is  an  accurate, 

convenient  and  fast  method  of  generating  hypothetical  data  sets 
for  the  activation  and  for  the  observed  photopeak  counts 
resulting  from  the  exposure  of  selected  nuclei  to  typical  flash 
and  cw  x-ray  spectra.  These  sets  of  data  formed  the  basis  for 
the  development  and  testing  of  the  second  program, 

UNFOLD,  a  deconvolution  procedure  capable  of  determining  the 
incident  x-ray  spectra  from  the  phoropeaks  observed  in  nuclear 
fluorescence  spectra  of  isomeric  decays  in  a  particular  set 
of  target  materials. 

ENFOLD 


ENFOLD  was  created  in  a  Lotus  1-2-3  spreadsheet.  This 
format  was  chosen  for  its  graphing,  editing  and  custom  menu 
features.  The  program  takes  into  account  many  factors  in  order 
to  accurately  simulate  data  measured  with  standard  scintillator 
spectrometer  techniques: 

1  .  Nuclear  -  Each  chosen  nucleus  possessing  an  isomeric 
state  has  its  own  set  of  nuclear  parameters,  such  as 
natural  abundance,  half-life  of  the  isomer,  dominant 
decay  fluorescence  energies,  and  relative  intensities 
of  the  signature  photons  .  Of  primary  importance  for 
each  nucleus  are  the  integrated  cross  sections  loD. 
and  energies  E^  of  the  gateways  through  which 
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photoactivation  proceeds  .  These  parameters  must  be 
manually  entered  by  the  operator. 

Detector  -  The  operator  has  an  option  for  each  of  the 

selected  nuclei  to  choose  with  which  of  3  typical 

detectors  the  sample  was  counted  in  the  simulation. 

The  options  are  between  a  3"  x  3"  solid  Nal(Tl) 

detector,  a  3"  x  3"  well  Nal(Tl)  detector,  or  a  2"  x  2" 

solid,  high  purity  germanium  detector.  The  choice  is 

made  through  a  custom  menu  in  the  spreadsheet  and 

provides  detector  efficiencies  at  the  fluorescence 
3 

energies  . 

Self -Absorption  -  Due  to  the  finite  thickness  of 
typical  samples,  some  self  absorption  of  fluorescence 
photons  occurs  and  a  suitable  correction  must  be  made. 
This  is  calculated  from  the  chosen  detector  and  typical 
sample  geometry  separately  by  Monte-Carlo  methods,  and 

4 

manually  entered  into  the  spreadsheet  . 

Decay  -  Correction  is  automatically  made  to  take  into 
account  the  decay  during  transport  of  samples  to  the 
detector,  counting  of  the  samples,  and  in  the  case  of 

4 

exposure  to  cw  x-ray  spectra,  during  the  irradiation  . 

Spectrum  -  The  operator  may  choose  to  simulate 

irradiation  of  the  samples  with  ■ a  variety  of  typical 
spectra,  or  may  enter  an  arbitrary  spectra.  Spectra 
included  in  the  spreadsheet  are  smoothed  transport  code 
evaluations  of  1.5  to  6  MeV  endpoint  cw  linac  spectra 
(normalized  to  the  output  of  a  Varian  Clinac  1800 
operating  at  300  Hz  in  a  6  MeV  photon  mode)5  and  of  6 
10  MeV  endpoint  flash  spectra  obtained  from  the 
DNA/ AURORA  nuclear  simulator5.  These  choices  are  made 
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through  custom  menus.  Once  a  spectrum  is  chosen,  the 
program  automatically  calculates  the  hypothetical, 
exact  activation  for  each  target  material  from 


r.:t  iaD.  ( 1  ) 

where  <?  ( )  is  tne  flux  a_  the  gateway  energy  . 
Spectra  may  be  plotted  in  several  different  formats 
through  a  custom  menu. 


6.  Random  deviates  -  Observed  photopeaks  always  contain 

statistical  variation  due  to  the  natural  Poisson 

distribution  of  the  decay  events.  The  operator  may 

choose  to  generate  observed  counts  directly  from  the 

exact  activation  or  to  include  the  statistical  effect 

by  generating  random  Poisson  deviates  about  the  exact 

7 

number  of  counts.  This  is  done  through  a  custom  menu. 


UNFOLD 

The  UNFOLD  program  was  written  in  MS-Fortran  in  order 
to  utilize  its  strong  debugging  capability.  The  deconvolution  of 
an  irradiating  spectrum  from  the  activations  present  in  target 
materials  depends  primarily  on  a  full  knowledge  of  the  cross 
sections  for  photoactivation  and  their  energy  locations  for  each 
material.  At  this  time,  this  information  is  not  available. 
However,  in  accordance  w^th  the  nuclear  systematic  study 
discussed  previously  it  is  possible  to  model  a  set  of  likely 
nuclear  parameters  by  the  following  assumptions: 

87  89 

1.  There  exist  a  set  of  four  isomeric  nuclei,  Sr,  Y, 
137  199 

Ba  and  Hg  which  are  activated  through  single 

gateways  located  between  1.5  and  3  MeV.  Within  this 
range,  the  gateways  were  placed  randomly  and  their 
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cross  sections  taken  directly  from  reference  2. 


There 
1  1  1 


exist 
1 13. 


a  set 
1 15. 


of  seven 
135 


77. 


isomeric  nuclei,  'Se, 

113  115_  135  1 79  ,  .  1 95_,  .  .  , 

Cd,  In,  In,  Ba,  Hf  and  Pt  which  are 

activated  through  two  gateways,  one  laying  between  1.5 

and  3  MeV,  and  the  other  laying  between  4  and  6  MeV. 

Within  each  range  the  gateways  were  located  randomly. 

For  each  material,  the  cross  section  of  the  upper 

gateway  was  taken  to  be  10  times  larger  than  that  of 

the  lower  gateway.  Values  for  the  cross  sections  were 

determined  from  reference  2  by  employing  eq.  1  and  the 

spectrum  of  the  device  used  in  these  experiments  to 

irradiate  the  samples,  the  aforementioned  linac 

operating  in  the  6  MeV  mode. 


The  modeling  procedure  was  to  determine  the  flux  at  each 
of  the  single  gateway  locations  and  then  to  fit  a  curve  to  these 
points,  using  general  least  squares  methods.  This  enabled  the 
flux  to  be  calculated  at  the  positions  of  the  lower  gateways  in 
each  multiple  gateway  material.  The  contribution  to  the  total 
activation  for  these  materials  from  the  lower  gateways  was  then 
be  removed  and  the  flux  at  the  location  of  the  upper  gateways 
found  using  eq.  1  . 


Three  examples  of  the  use  of  UNFOLD  are  included  here  to 
illustrate  the  degree  to  which  this  approach  has  been  successful 
in  deconvolving  the  irradiating  spectrum.  Example  1  shows  the 
results  of  a  3  MeV  endpoint  irradiating  linac  spectrum.  In  this 
example  none  of  the  upper  gateways  of  the  multiple  gateway 
materials  were  accessed.  Example  2  shows  the  results  of  a  4.5 
MeV  end  point  linac  spectrum  in  which  some,  but  not  all,  of  the 
upper  gateways  were  accessed.  Finally,  example  3  displays  the 
unfolding  of  a  6  MeV  end  point  linac  spectrum  in  which  all 
gateways  are  accessed.  In  each  case  the  extent  to  which  the 
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irradiating  x-ray  spectrum  was  reproduced  at  the  gateway 
locations  was  remarkable.  The  high  degree  of  success  achieved 
using  this  model  is  indicative  of  the  potential  of  the  XAN 
technique.  The  problem  remains  to  experimentally  determine  the 
locations  and  cross  sections  for  photoactivation  for  a  wide  range 
of  nuclei,  a  problem  which  cannot  be  adequately  addressed  at  this 
time . 
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EXAMFLE  1 


MeV 


EE AT I  OK  PROGRAM  -  (ENFOLD)  -  Produces  observed  counts  represent- 
activationc  present  in  hypothetical  samples  due  to  e  particular 
ei.  input  >:-ray  spectrum 

Genera.  Instructions 

ice:  c f  input  spectrum,  detector  efficiencies,  type  or  random 
deviates,  and  spectral  plot  (if  desired'  are  made  through  menus. 

Lpectrun  menu  -  L inac , Arbi trar y  or  Aurora  -  Alt -D 
Detector  menu  -  Nal  well  or  solid  or  HPGe  well  -  7iit-D 


r.  d  o  r.  a  e  v  1  a  t  e  s 

me 

nu  - 

Exact  or 

Poisson  -  Alt-F. 

ectra-  piot  menu 

inear f loo 

or  MeV/MeV  -  Alt 

Linac  EEnd ! M 

eV  i 

= 

3 

Isomer  1 :  3"  x 

*• 

Nal 

well 

I some r  2  :  ?"  x 

»> 

Nal 

well 

Isomer',:  3''  x 

*>  M 

Nal 

well 

Isomer  4:  3"  x 

“s  fi 

Nal 

well 

’ :  P :  l  s cor. 


Spectral  Distribution 


Energy [ Me 


Phi [ i/cm" 2-keV-sec ] 


G.  0  6 : 
0.196 
C  .  313 
0 .436 
0.563 
0.686 
0.813 
0.336 
1.063 
1.18  8 

1.313 
1.438 

1.563 
1  .688 

1 .813 

1.938 
1.062 
2 . 188 

2 .313 
2  .438 

2.563 
2 .688 

2 .813 

2.938 


1 . 47E+06 
4 . 79E+09 
1 . 69E+10 
1 . 67E+10 
1 . 64E+10 
1 . 55E+10 
1 . 38E+10 
1 . 06E+10 
8 . 69E+09 
7 . 23E+09 
6 . 08E+09 
5.16E+09 
4 . 41E+09 ' 

3 . 80E+09 
3 . 28E+09 
2 . 86E+09 
2 . 49E+09 
2.10E+09 
1 . 77E+09 
1 . 42E+09 
1 . 09E+09 

8 . 02E+08  48 

5.63E+08 

3.77E+08 


Target  Isomer  Information 
Single  gateway  nuclei 


Isomer  1 

8  7Sr 

Isomer  2 

8  91 

NTar gets 

=  1.0E+20 

Ntar gets 

=  1.0E+22 

EGatelheV]  SG[ 

1 C  ~ - 2  9  cm  ~  2 - k  e  V ] 

EGate [MeV] 

SG  [  10  '  -  29cn," 

2.6c 

20  0  8 

2 . 79 

59  2 

0 

0 

0 

0 

0 

0 

0 

0 

f 

0 

0 

0 

c 

n 

0 

0 

0 

0 

0 

0 

Exact Act [ 1/sec  ]  = 

1 . 73E-17 

ExactAct [ 1/ sec  )  = 

3 . 56E-18 

Tl/2 [ s ;  = 

1 . 01E+04 

Tl/2 [ s  ]  = 

1.61E+01 

IrrTis]  = 

300 

IrrT[s]  = 

300 

TTls] 

60 

TTls] 

60 

CntTls]  = 

300 

CntTls]  = 

300 

DecCorri= 

2 . 03E-02 

DecCorrl= 

7 . 55E-02 

DecCor r  2  - 

9 .90E-01 

DecCorr 2= 

7 . 74E-02 

DeadTCor - 

1 . OOE+OO 

DeadTCor = 

1 .00E+00 

E 1  [  k  e  V  ]  = 

388 . 4 

EllkeV]  = 

909.15 

Intens  = 

0 .823 

Intens  = 

0 . 9914 

SelfAbs  = 

1.000 

SelfAbs  = 

1.000 

DetEff 

0 .339 

DetEff  = 

0.102 

ODsCtsl  = 

2937 

■ObsCtsl  = 

6255 

StandDev- 

54 

StandDev= 

79 

E  2  [  k  e  V  ]  = 

0 

E  2  [  k  e  V }  = 

0 

I  n  t  n  s  = 

0 

Intens  = 

0 

SelfAoc  =• 

0 .000 

SelfAbs  = 

0.000 

DetEff  = 

0 .000 

DetEff  = 

0.000 

ObsCts2  = 

0 

ObsCts2  = 

0 

StaridDev  = 

0 

StandDev= 

0 
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Isomer  3-  13  7Ba 

NTar gets  -  i.OE+19 

EGa  te  l  MeV  3  SGllO~-29cm"2-keV] 

I  .  ?  I  2203 

U  0 

ri  0 

0  0 

0  0 

r  n 


ExactAct [i/ssc  i  = 

6 . 01E-17 

Tl/2 I s :  = 

I . 53E+01 

IrrTls ]  - 

30  0 

TTls! 

6  0 

CntTls]  = 

3  0  0 

DecCor r 1= 

5 . 66E-01 

DecCorr  2- 

5 . 47E-01 

DeadTCor  = 

1 . 00E+00 

EllkeV;  - 

661.66 

Intens 

0 .901 

SelfAbs  = 

1.000 

DetSfl  = 

G  .  181 

ObsCtsi  = 

9089 

StandDev- 

95 

E2 [ keV ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0.000 

DetEff  = 

C  .000 

OosCtsl  = 

0 

StandDe v= 

c 

Isomer-;-  199Hg 
Ntargets=  I.OE+19 

EGa  te [ Me V ]  SG [ 10' -29cnT  2 -keV ] 


2 .83 

3698 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Exact Ac 1 1 1 /sec  )  = 

1 . 9  7E-1'7 

Tl/2 i s  3  = 

2 . 5  6E  +  0  3 

IrrTls]  = 

300 

TTls] 

60 

CntT l s ’  = 

300 

DecCorr 1= 

7 . 69E-02 

DecCorr  2  = 

9  .  60E-01 

DeadTCor = 

1 . OOE+OO 

EllkeV]  = 

156 . 376 

Intens  - 

0 . 53 

SelfAbs  = 

1.000 

DetEff  = 

0 . 577 

ObsCtsi  = 

1372 

StandDev= 

37 

E2 [ keV ]  = 

0 

Intens  = 

0 

SelfAbs  - 

0 . 000 

DetEff  = 

0 .000 

ObsCts2  = 

0 

StandDev= 

U 
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Target  Isomer  Information 
Multiple  gateway  nuclei 


Isomer  1 

-  7  7Se 

1 somer  2 

-  lllCd 

NTargets 

=  1.0K421 

Ntar gets 

=  1.0E+20 

E  G  a  t  e  [  M  e  V  )  SG[lC'-29cm~2-keV] 

EGate i MsV 1 

S  G [ 1C  ~ -  29cm" 

1.75 

289 

1 .  G  2 

53 

.  - 
*1  .  X  .1 

2991 

4.26 

529 

0 

0 

0 

0 

c 

C 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

ExactAct[l/sec]- 

1 . 02E-17 

ExactAct[l/sec]= 

2 . 18E-18 

T 1  /  2  [  s  )  = 

1 . 75E  +  01 

Tl/2 l s ]  = 

2 . 92E  +  03 

lrrT[s)  = 

300 

IrrTls]  = 

300 

TTts] 

6  C 

TT  [  s  i 

60 

CntTIs;  =• 

30  0 

CntTIs)  = 

300 

DecCcrrl= 

3  .  °2E-02 

DecCorr 1= 

6  .  79E-02 

DecCorr  2  = 

8 . 39E- u  2 

DecCorr  2  = 

9 . 65E-01 

DeadTCor  = 

1 . OOE+OC 

DeadTCor = 

1 . C0E+00 

El [ KeV ]  = 

1  o  1  .32 

El(keV)  = 

245.384 

Intens  = 

G  .  524 

Intens  = 

0  =<4 

SelfAbs  = 

1 . 000 

SelfAbs  = 

1.000 

DetEff  = 

0  .  572 

DetEff  = 

0.472 

ObsCtsl  = 

7171 

ObsCtsl  = 

1871 

StandDev= 

85 

StandDev= 

43 

E2[KeV]  = 

G 

E2 l keV }  - 

0 

Intns  = 

C 

Intens  = 

0 

SelfAbs  = 

C  .  0  0  0 

SelfAbs  - 

0 . 000 

DetEff  = 

C  .  0  0  0 

DetEff  = 

0 . 000 

ObsCtsl  = 

c 

ObsCtsl  = 

0 

StandDev= 

0 

StandDev= 

0 

-keV  ] 
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Isomerl-  1131n 
NTarqets=  5.0E+2C 


Isomer4-  1151n 
Ntargets=  1.0E+21 


EGate . MtV ]  S  G  [ 

1C  '  - 29cm~2-keV ] 

EGate [ MeV ] 

SG[10"-29crri" 

1.96 

80 

2  .  46 

123 

4  .  51 

79  7 

4.82 

1228 

0 

C 

0 

0 

r 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Exact Act [ 1/sec  i  = 

2  .  2  3  E  - 1 6 

Exact Act [ 1/sec  1  = 

1  .  67E-18 

Ti / 2 i s  ]  - 

•  .  97E+03 

Tl/2 [ s  1  = 

1 . 57E+04 

IrrTIs]  = 

3  0  C 

IrrTIs]  = 

300 

TT  l  s  ] 

60 

TTls] 

60 

CntTls]  = 

300 

CntTls]  = 

300 

DecCorr 1= 

3  .  40E-02 

DecCorr 1= 

1 . 32E-02 

DecCorr 2 = 

5  .  83E-01 

DecCorr  2  = 

9 . 9  3E- 0 1 

DeadTCor - 

1 . 00E+00 

DeadTCor = 

1 . 00E+00 

EllkeV]  = 

391.69 

EllkeV]  - 

336.26 

Intent 

0.642 

In tens  = 

0 .458 

SelfAos  = 

1.000 

SelfAbs  = 

1 .000 

DetEff 

0 .337 

DetEff  = 

0.383 

Onsets  1  = 

2  4  8  6 

ObsCtsl  = 

116  5 

StandDev= 

50 

StandDev= 

34 

E  2  [  k  e  V  ]  = 

0 

E  2  [  k  e  V  ]  = 

0 

Intens  = 

0 

Intens  = 

0 

SelfAos  = 

0 .000 

SelfAbs  = 

0 . 000 

De tEf  f  = 

0.000 

DetEff 

0.000 

OocCtsl  = 

0 

ObsCtsr  = 

0 

StandDev= 

0 

StandDev= 

0 
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Target  Isomer  Information 
Multiple  gateway  nuclei 


Isomer  5-  125&s  Isomer6-  173Hf 

Mar  gets=  1.0E-t22  Ntargets=  1.0E+20 

EGatelMeV;  SG ! 1 2 '  -23cm" 2 -keV ]  EGatetMeV]  SG [ 1 C ' - 29 cm" 2 -ke V ; 


n  c 

211 

2.08 

494 

4.21 

2110 

4.69 

49  4  2 

c 

0 

0 

0 

l 

C 

0 

0 

c, 

v' 

0 

0 

0 

0 

0 

0 

Exact Act [ 1 / se c  j  = 

7  .  62E-19 

Exact Act [ 1 /sec  j  = 

1 . 21E-17 

71/2  is]  = 

1 . 03E+05 

Tl/2 ( s  )  = 

1 . 87E+01 

1  r  r  T  i  s  ]  = 

300 

IrrTi s ]  = 

300 

TTls] 

60 

TT  [  s  ] 

60 

C  n  t  T  t  s  ]  = 

300 

CntTls]  = 

300 

DecCor  r  1  - 

2 . 01E-03 

DecCorr 1= 

1 . 08E-01 

DecCorr 2= 

5 . 99E-01 

DecCorr  2  = 

8 . 98E-02 

DeacTCcr = 

1 . 00E+0O 

DeadTCor  = 

1 .00E+00 

e: ; ke v  ]  = 

266 . 27 

ElikeV]  = 

214 . 31 

Intens  = 

0.156 

Intens  = 

0.952 

SelfAbs  - 

1.000 

SelfAbs  = 

1.000 

DetEff 

C  .  448 

DetEff 

0 . 507 

ODtCtcl  = 

306 

ObsCtsl  = 

1713 

StandDev= 

17 

S tandDev= 

41 

E  2  [  k  e  V  ;  - 

0 

E  2  [  k  e  V  j  - 

0 

lntns  = 

0 

Intens  = 

0 

SelfAcs  = 

0 . 0  0  0 

SelfAbs  = 

0 .000 

DetEff 

0.000 

DetEff  = 

C  .000 

ObsCtsl  = 

0 

ObsCts2  = 

0 

StandDe v= 

0 

StandDev= 

0 
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1 s  one :7- 

19  5Pt 

NTa raers- 

3- .  0  E  +  2  2 

Elate 

; Mev ;  sg; ic 

~  -  29 crrt “ 

- 

2  22 
^  ^  ^  — 

.  *1 

L' 

^  J 

0 

r\ 

c 

0 

(' 

0 

Exact 

Art! 1/ sec;- 

. 26E-1E 

Tl/2  [ 

c  ]  —  ^ 

.  47E+05 

IrrT  l 

S  j  ~ 

300 

TT  i  s  ] 

= 

6  0 

CntTl 

s  l  — 

3  0  0 

DecCorrl=  5 

.  98E-04 

DecCcrr2=  1 

.  00E  +  00 

DeadTCor=  i 

. 00E+00 

E 1  [  k  e  V  J  = 

1 2  S  .  7  6 

1  nten 

C-  IT 

0.02835 

SelfAbs  = 

1  .000 

DetEf 

t-  - 

G  .616 

ObsCt 

si  — 

109  5 

Stand 

Lev- 

3  3 

E2 [ xeV  ;  =  0 
Inrens  =  C 
SelfAbs  =  O.OOC 
DetEfr  =  O.OOC* 

Ods2cs2=  0 
StandDev=  0 


C:  \ NEWPORT! PR OGR AMS  >unf  old 
Enter  1  for  a  cv  exposure  : 

Enter  the  end  point  energy  [MeV]  of  the  x-ray  device: 

Enter  Tlrr,  TTrav  and  TCnt  [seel  for  the  87Sr  sample: 

3  0  C  60  30 C 

Enter  the  numDer  of  target  nuclei  for  the  87Sr  sample: 

1E2C 

Enter  total  observed  and  bkgnd  counts  in  the  386  keV  peak: 
29  2  ”  0 

Enter  C,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  tne  se 1 f -absor pt i on : 

Enter  Tlrr,  TTrav  and  TCnt  iseci  for  the  89Y  sample: 

3  0  C  60  300 

Enter  the  number  of  target  nuclei  for  the  89Y  sample: 

1E2  2 

Enter  total  observed  and  bkqnd  counts  in  the  909  keV  peak: 
•32  5  5  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  sel f -absor pt i on : 


Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  137Ba  sample: 

3  0  C  60  300 

Enter  the  number  of  target  nuclei  for  the  137Ba  sample: 

1E19 

Enter  total  observed  and  bkgnd  counts  in  the  661  keV  peak: 

9089  0 

Enter  0,  1  or  2  for  a  Nal  solid,  veil,  or  HPGe  detector: 

Enter  the  se 1 f -absor pt i on : 

JL 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  199Hg  sample: 

300  60  300 

Enter  the  number  of  target  nuclei  for  the  199Hg  sample: 

1E1 9 

Enter  coral  observec  and  bkgnd  counts  in  the  156  keV  peak: 

1372  0 

Enter  0,  1  or  2  for  a  Nal  solid,  veil,  or  HPGe  detector: 

Enter  the  self-absorption: 

Enter  1  to  examine  the  fit  employed  in  the  single  gatevay  region: 

The  coefficients  of  the  pover  series  are: 

r.  =  1,  A  =  5.068  59  5E  +  09 

r.  =  2,  A  =  2 . 81830 2E  +  09 

r.  =  3,  A  =  -3.08968  2E  +  09 

r.  =  4,  A  =  5 . 397999E+08 

Chi  squared  =  1.512296 


Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  77Se  sample: 
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IOC 

6C 

300 

Ent 

the  numDer  of  target 

nuclei 

for 

the 

77Se 

sample : 

Ent 

7^7 

e.  '• 

1  c 

—  Z 

total  oDse.rved  and  bkgnd  counts 

in 

the  162  keV  peak 

Ent 

C ,  _l  or  2  for  a  Nal 

solid. 

well 

,  or 

HPGe 

detector : 

Ent 

6*  Z 

the  se 1 r -absor pt i on : 

Ent 

Tlrr,  TTrav  and  TCnt 

[  sec  ] 

for 

the 

11  lCd 

sample : 

3  o  r 

6  0 

30  u 

^  r*  t. 

d  >* 

the  number  of  target 

nuclei 

for 

the 

lllCd 

sample : 

TTf' 


Enter  total  observed  and  bkgnd  counts  in  the  245  keV  peak 
1 E  ~  1  C 

Enter  C,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  self -absorpti on : 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  113ln  sample: 

3  0  C  60  300 

Enter  the  numDer  of  target  nuclei  for  the  113In  sample: 
5E2  0 

Enter  total  observed  and  bkgnd  counts  in  the  392  keV  peak 

2  4  6  c  0 

Enter  C,  I  cr  2  for  a  Nal  solid,  well,  or  HPGe  detector: 


Enter  the  self -absorption : 


Enter 

2  0  C  6  ( 

Tlrr,  TTrav  and  TCnt 
i  30C 

[  sec  ] 

for 

the 

1 1 5 1  n 

sample  : 

Enter 
1E2  1 

tne  number  of  target 

nuclei 

for 

the 

1 1 5 1  n 

sample : 

Enter  total  observed  and  bkgnd  counts  in  the  336  keV  peak 
116  5  0 

Enter  C,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  self-absorpti on : 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  135Ba  sample: 

3 0 C  60  300 

Enter  the  number  of  taraet.  nuclei  for  the  135Ba  sample: 
1E2  2 

Enter  total  observed  and  bkgnd  counts  in  the  268  keV  peak 
30  6  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  self-absorption: 


(rt  (_*  t1’ 


Enter  TIrr,  TTrav  ano  TCnt  [sec]  for  the  179Hf  sample: 

3  0  :  61  3  0  c 

r.ter  me  nuir.De:  of  target  nuclei  for  the  179Hf  sample: 

E  <.  3 

nter  rota]  observed  and  bkgnd  counts  in  the  214  keV  peak: 
171 :  C 

Enter  C >,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector:- 
Enter  the  self-absorption: 

Enter  TIrr,  TTrav  and  TCnt  isec]  for  the  195Pt  sample: 

33:  60  300 

Enter  tne  number  of  target  nuclei  for  the  195Pt  sample: 
5E21 

Enter  total  observed  and  bxqnd  counts  in  the  130  keV  peak: 
109:  0 

Enter  0 ,  *  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  se 1 r -abs cr pt i on : 


Tnt  values  of  gateway  energy,  flux  and  standard  deviaton  are: 


2 . 660000 

8 . 7C1814E  +  08 

1. 605676E+07 

2.790000 

5.99  3  56 1E+  0  8 

7578291.000000 

1 .980000 

2.725342E+09 

2 . 858663E+07 

2.830000 

5 . 480612E+08 

1. 479626E+07 

1 .750000 

3 . 56499 0E+09 

4 . 209868E+07 

1 .620000 

4 . 05039  2E  +  09 

9 . 363972E+07 

1.960000 

2 . 866556E+09 

5.74923  3E  +  0  7 

2 .460000 

1 . 37  7  69  7E+09 

4 . 036368E+07 

2 .950000 

3 . 447520E+08 

1 . 970816E+07 

2 .080000 

2 . 470555E+09 

5 . 969197E+07 

1 .820000 

3 . 07899 8E+09 

9 . 304698E+07 

2:  \ NEWF OR T\ PROGRAMS  > 
C : \ NEWPORT \ PROGRAMS  > 
C : \ NEWFORTXPROGRAMS  > 
C : \NEWFORT\PROGRAMS> 
C : \NEWFORT\PROGRAMS > 
C : \ NEWFORTXPROGRAMS > 
C : \ NEWFORTXPROGRAMS  > 
C : \ NEWFORTXPROGRAMS  > 
C : \ NEWFORTXPROGRAMS  > 
S : \NEWFORT\PROGRAMS> 
C : \ NEWFORTXPROGRAMS  > 


57 


Fit  b'i 


EXAMPLE 


n 


DATA  GENERATION  PROGRAM  -  (ENFOLD)  -  Produces  observed  counts  represent¬ 
ing  activations  present  in  hypothetical  samples  due  to  a  particular 
give;,  input  x-ray  spectrum 


General  Instructions 

Cm  ices  of  input  spectrum.,  detector  efficiencies,  type  of  random 

deviates,  and  spectral  plot  (if  desired)  are  made  tnrouqh  menus. 

Spectrum,  menu  -  Linac,  ArDitrary  or  Aurora  -  Alt-S 
Detector  menu  -  Nal  well  or  solid  or  HPGe  well  -  Alt-D 
Random  deviates  menu  -  Exact  or  Poisson  -  Alt-R 
Spectral  plot  menu  -  Linear, Log  or  MeV/MeV  -  Alt-P 


£  peter  urn. : 
Detector  : 

Rant  uev: 


j_,mac 

1 somer  1 
Isomer  2 
1 s omer  3 
Isomer  4 

Po  i  ss  or. 


EEnd [ MeV  j  = 


3  " 

X 

3  " 

3  " 

X 

?  M 

?  M 

X 

r> »» 

o  M 

V 

v 

4 . 5 

Nal  well 
Nal  well 
Nal  well 
Nal  well 


Spectral  Distribution 


Energy [ MeV 1 


Phi  [  1/cm" 2-k.eV-sec  ] 


C  .  094 

6 . 56E+07 

0 .281 

2 . 14E+09 

0.469 

7 . 52E+09 

0.6  56 

7 . 45E+09 

0.644 

7 . 33E+09 

_  .  C  2  : 

6 . 91E+09 

1 . 2  1  Ci 

6 . 15E+09 

1.406 

4 . 72E+09 

1.594 

3  .  87E  +  09 

1.781 

3 . 22E+09 

1.969 

2 . 71E+09 

2.156 

2 . 30E+09 

2.344 

1 .97E+09 

2 . 531 

1 . 69E+09 

2 .719 

1 . 46E+09 

2.906 

1 . 27E+09 

3.094 

1 . 11E+09 

3 .281 

9 . 36E+08 

3.469 

7  .  87E  +  08 

3.656 

6  .  34E+08 

3.844 

4  .  87E  +  08 

4.031 

3 . 58E+08 

4 . 219 

2 . 51E+08 

4.406 

1 . 68E+08 
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Target  Isomer  Information 
Single  gateway  nuclei 


Isomer  1-  87Sr  lsomer2-  89Y 

KTar  gets=  ^.0E+20  Ntargets=  1.0E  +  22 


EGa'-fc  [  MeV  ; 


SG  [  1 C  '  -  2  9  err. '  2  -  keV  ]  EGatetMeV]  SG  [  1  0  ~  -  29cm"  2 


_  .  6  6 

2  0  0  6 

0 

0 

0 

0 

0 

0 

u 

0 

0 

0 

Exact Ac 

t [ 1/sec i = 

2  .  07E-17 

Tl/2ls] 

= 

1 . 01E+04 

I rrTl s  ] 

= 

300 

TT  i  s  ] 

= 

60 

CntT [ s  ] 

300 

DecCor r 

1  = 

2 . 03E-02 

DecCorr 

2  = 

9 . 90E-01 

DeadTCor = 

1 . 00E+00 

EltkeV] 

386.4 

1 ntens 

= 

0 .023 

Se 1 f Abs 

= 

1 .000 

DetEf f 

- 

0 .339 

ObsCts 1 

5235 

StandDev= 

72 

E  2  t  k  e  V  ]  = 

0 

lntnr  - 

0 

SelfAos  = 

0 .000 

DetEff  = 

0 . 000 

ObsCtsr  = 

0 

StandDev= 

0 

2 . 79 

592 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ExactAct[l/sec]= 

8  .  21E-18 

Tl/2  l  s  ]  = 

1 . 61E+01 

IrrTls]  = 

300 

TT  [  s  ] 

60 

CntTls]  = 

300 

DecCorr 1= 

7 . 55E-02 

DecCor r  2  = 

7 . 74E-02 

DeadTCor = 

1 .00E+00 

El [ keV ]  = 

909 .15 

Intens  = 

0.9914 

SeifAbs  = 

1 .000 

DetEff  = 

0 .102 

ObsCtsl  = 

14590 

StandDev= 

121 

E2 [ keV ]  = 

0 

Intens  = 

0 

SeifAbs  = 

0 .000 

DetEff  = 

0.000 

ObsCts2  = 

0 

StandDev= 

0 

-keV  ] 


60 


Isomer3-  137Ba 
NTarqets-  i.OE+19 

EGacelMeV;  SG110'-29cm~2-keV] 


1.91 

2203 

0 

0 

0 

r 

0 

r. 

L 

0 

U 

0 

r 

0 

Exact Act [ 1/sec  j  = 

5.91E-17 

T 1  /  2  i  s  1  = 

1 . 53E+02 

IrrT [ s ;  = 

300 

TTls! 

60 

CntT [ s ]  - 

300 

DecCor r 1= 

5 . 66E-01 

DecCor r  2  = 

5 . 47E-01 

BeadTCor- 

1 . 00E+00 

EllkeV]  = 

661.66 

Intens  = 

0.901 

SelrAbs  = 

1 . 000 

DetEff  = 

0.181 

ObsCtsl  = 

87  4  4 

StandDe v= 

94 

E  2  [  k  e  V  ]  = 

0 

Intens  = 

0 

SelfADS  = 

0.000 

DetEff  = 

0 . 000 

ObsCtsl  = 

0 

StandDev- 

0 

Isomer4-  199Hg 
Ntargets=  I.OE+19 

EGa  te [ Me V ]  SG 1 1 0 ~ - 29cm~ 2 -ke V ] 


2.83 

3698 

0 

0 

0 

0 

0 

C 

0 

0 

0 

0 

ExactAct [ 1/sec  ]  - 

4 . 98E-17 

Ti/2 l s ]  = 

2 . 56E+03 

IrrTIs]  = 

300 

TTls  3 

60 

CntT [ s  ]  = 

300 

DecCorr 1= 

7 . 69E-02 

DecCor r  2- 

9 . 60E-01 

DeadTCor  = 

1 . 00E+00 

ElfkeV]  = 

158 . 376 

Intens  = 

0.53 

SelfAbs  = 

1 .000 

DetEff  = 

0 . 577 

ObsCtsl  - 

3447 

StandDev= 

59 

E2 [ keV ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0 .000 

DetEff  = 

0 . 000 

ObsCts2  = 

0 

StandDe v= 

0 
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Target  Isomer  Information 
Multiple  gateway  nuclei 


Isomeri-  77Se  Isomer2-  lllCd 

NTargets-  1.0E+21  Ntargets=  1.0E+20 

EGatetMeV]  SG  [  1 0  "  -  29cm''  2 -keV  ]  EGatetMeV]  SG  [  10 '  -  29cm  ~  2 


1  . 

289 

A  *  "S 

289  1 

r : 

0 

0 

0 

0 

G 

0 

Exact Act [ ‘ /sec j  = 

1 . 84E-17 

Tl/2ls]  = 

1 . 75E+01 

IrrTts]  = 

300 

TT  [  s  ] 

60 

CntT [ s  j  - 

2  n  r: 

DecCor r 1 = 

S . 22E-02 

DecCor r  2  = 

6 . 39E-02 

DeadTCor  = 

1 . 00E+00 

EllkeV]  = 

161 .92 

I  n  *2  o  r]  c  — 

0 . 524 

Self Ads  = 

1 . 000 

De  tEf  £  = 

0 . 572 

OosCtsl  = 

12838 

StandDev= 

113 

E  2  [  k  e  7  ;  = 

0 

Intns  = 

0 

SelfAbs  = 

0.000 

DetEff 

0 . 000 

ObsCtsl  - 

0 

StandDe v= 

0 

1 .62 

53 

4 . 26 

529 

0 

0 

0 

0 

0 

0 

0 

0 

Exact Act [ 1/sec ] = 

3.22E-18 

Tl/2[s]  = 

2 . 9  2E  +  0  3 

IrrTts]  = 

300 

TTts] 

60 

CntT [ s  3  = 

300 

DecCorr 1= 

6 . 79E-02 

DecCorr 2= 

9 . 65E-01 

DeadTCor = 

1 .00E+00 

EllkeV]  = 

245.384 

Intens  = 

0.94 

SelfAbs  = 

1.000 

DetEff  = 

0 .472 

ObsCtsl  = 

2760 

StandDev= 

53 

E  2  [  k  e  V  ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0 .000 

DetEff  = 

0.000 

ObsCts2  = 

0 

StandDev= 

0 

-keV  ] 
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Isomer  3-  11 3 In  Isomer  4-  1151 n 

KTargets=  5.0E+2G  Ntargets=  1.0E+21 

EGate  [  MeV  3  SG  :  1 0 29cnT  2 -ke  V  ]  EGatelMeV]  SGI  10~-29cni~2-keV  j 


1.46 

80 

2.46 

123 

A  Z  “ 

M  .  J. 

797 

4.82 

1228 

C 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Exact Act [ 1/sec  j  = 

2 . 18E-18 

ExactActIl/sec]= 

2 . 20E-18 

Ti/2 [ s  3  = 

5 . 97E+03 

Tl/2 [ s )  = 

1 . 57E+04 

IrrTis]  = 

300 

IrrTis]  = 

300 

TTls] 

60 

TT  [  s  ] 

60 

CntTlsJ  = 

300 

CntTIs]  = 

300 

DecCorr 1= 

3 . 40E-02 

DecCor r 1= 

1. 32E-02 

DecCorr  2  = 

9 . 83E-01 

DecCor  r  2  = 

9 . 9  3E-01 

DeadTCcr = 

1.00E+00 

DeadTCor = 

1 . 00E+00 

El [ keV  3  = 

391.69 

EllkeV]  = 

336.26 

Intens  = 

0 .642 

Intens  = 

0 .458 

SelfAbs  = 

1 . 000 

SelfAbs  = 

1.000 

DetEff  = 

0.337 

DetEff  = 

0.383 

OcsCtsl  = 

2404 

ObsCtsl  = 

1474 

StandDe v  = 

49 

StandDev= 

38 

E2 [ keV ]  = 

0 

E2 I keV  ]  = 

0 

Intens  = 

0 

Intens  = 

0 

SelfAos  = 

0 .000 

SelfAbs  = 

0 . 000 

DetEff  = 

0. 000 

DetEff  = 

0.000 

ObsCts2  = 

0 

ObsCts2  = 

0 

StandDev= 

0 

StandDev= 

0 
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Target  Isomer  Information 
Multiple  gateway  nuclei 


Isomer  5-  12  5Ba 

NTarqets=  1.0E+22 


EGate ! MeV ] 

SG[ lC~-29cnr 

2  .  9  r- 

211 

4.21 

2110 

0 

0 

C 

0 

r 

0 

C' 

r> 

o 

ExactAct [ 1/se 

cJ=  7.99E-18 

Tl/2 [ s ]  = 

1 . 03E  +  05 

IrrTls]  = 

300 

TTls] 

60 

CntTls]  = 

300 

DecCor r 1= 

2 . 01E-03 

DecCorr  2  = 

9 . 99E-01 

DeadTCor  = 

1 . 00E+00 

EllkeV]  = 

268 . 27 

Intens  = 

C  .156 

SelfAbs  = 

1 .000 

DetEff  = 

0.448 

ObsCtsl  - 

3  3  57 

StandDev= 

58 

E2 [ ke V  j  = 

0 

1  n  c  n  s  - 

0 

SelfAbs  = 

C  .  000 

DetEff 

0  .  ooc 

OosCts2  = 

0 

S  tandue v  = 

0 

Isomer6-  179Hf 
Ntargets=  1.0E-I20 

EGate  [  MeV  j  SG  [  1 C  ~  - 29cm~  2 -ke V  ] 


2.08 

494 

4  .  69 

4942 

0 

0 

0 

0 

0 

0 

0 

0 

ExactAct[l/sec]= 

1.21E-17 

Tl/2 [ s ]  = 

1 . 87E+01 

IrrTls]  = 

300 

TTls] 

60 

CntTls]  = 

300 

DecCorr 1= 

1 . 08E-01 

DecCorr  2  = 

8 . 9  8E- 0  2 

DeadTCor  = 

1 . OOE+OO 

EllkeV]  = 

214.31 

Intens  = 

0.952 

SelfAbs  = 

1 .000 

DetEff  = 

0 . 507 

ObsCtsl  = 

1752 

StandDe v= 

42 

E2 l keV ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0 .000 

DetEff  = 

0 .000 

ObsCts  2  = 

0 

StandDev= 

0 
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Isomer?-  195Pt 
NTargets-  5.0E+22 

EGatelMeVj  3G  [  1 0  29cm''  2-keV  ] 

223 
2  3  2 
0 
0 
0 
0 

Exactiioo  [  i/sec  i  =  G.92E-18 

3 . 4  7E  +  G5 
300 
60 
300 
5 . 98E-04 
1 . 00E+00 
1 . 00E+00 


El ! ke V ]  = 

129  .76 

|  rj-pnc  - 

C .02885 

SelfAbs  = 

1 .000 

DetEff  = 

C  .616 

CDSCtSl  - 

1099 

SrandDe v= 

33 

E  2  [  k  e  V  i  = 

0 

Intens  = 

0 

SelfAbs  = 

0.000 

LetEff 

0  .  000 

ODsOtsl  = 

0 

StandDev= 

0 

TT  i  s  ]  = 

CntT [ s  ]  = 
DecCorr  1  = 
DecGor  r  2  = 
DeadTGor = 


2  .  G  2 
.  4 1 
C 
0 
0 
0 


:■ :  \NEWFORT  \  PROGRAMS  >  unfold 
Enter  1  ter  6  tv  exposure: 

Enter  tne  enc  point  energy  [MeVJ  of  the  x-ray  device: 

A  ,  tx 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  87Sr  sample: 

3  0  C  6  C  3  0  C 

Enter  the  number  or  target  nuclei  for  the  87Sr  sample: 

1E2C 

Enter  total  observed  and  bkanc  counts  in  the  388  keV  pea);: 

5  :  3  f  c 

Enter  C,  1  or  1  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  tne  se 1 : -absorpt i on : 

Enter  Tlrr,  TTrav  and  TCnt  [seci  for  the  89Y  sample: 

3  0  C  G  C  3  0  C 

Enter  the  number  of  tarqet  nuclei  for  the  89Y  sample: 
iE2: 

Enter  total  observed  and  bkanc  counts  in  the  909  keV  peak: 
14590  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  self -absorpti on: 


Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  137Ba  sample: 

300  60  300 

Enter  the  number  of  target  nuclei  for  the  137Ba  sample: 

IE  19 

Enter  total  observed  and  bkand  counts  in  the  661  keV  peak: 
3  7  4  4  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 

■» 

i 

Enter  the  self -absorpt i on : 

Enter  Tlrr,  TTrav  anc  TCnt  [sec]  for  the  199Hg  sample: 

3  0  C  60  300 

Enter  tne  number  of  target  nuclei  for  the  199Hg  sample: 

1E1 c 

Enter  total  observed  and  bkgnc  counts  in  the  158  keV  peak: 
3  4  4^  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 

± 

Enter  the  sel f -absorpti on : 


Enter  1  to  examine  the  fit  employed  in  the  single  gateway  region 


Tne  coefficients  of  the  power  series  are: 
n  =  1,  A  =  3 . 657987E+09 

n  =  2,  A  =  2.11557  4E+  09 

n  =  3,  A  =  -2 . 028464E+09 

n  =  4,  A  =  3 . 514063E+08 

Chi  squared  =  4.894636E-01 


Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  77Se  sample: 


3  0  C  6  C 

3  0  0 

Enter 
1E2  : 

tne 

number  of  target 

nuclei 

for 

the 

7  7Se 

sample  : 

Enter 

t  eta 

1  observed  and  bkgnd  counts 

i  n 

the  16 

2  keV  pea 

12  6  3  8 

r 

Enter 

r  ■ 

/  - 

or  2  for  a  Nal 

solid. 

well 

,  or 

HPGe 

detector  : 

Enter 

tn  l 

self-absorption: 

V  n  t_  6: 1 

,  TTrav  and  TCnt 

[  sec  ] 

f  or 

the 

lllCd 

sample • 

3  0  0  o  C 

5  0  C 

Enter 

tne 

number  of  target 

nuclei 

for 

the 

illCd 

sample  : 

IE 2  C 

Enter  total  observed  and  bkgnd  counts  in  the  24:.  keV  pear. 

2  7  u  0  C 

Enter  2,  1  cr  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  tne  se 1 f -abs or pt i on : 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  113Ir:  sample: 

3  0  0  6  0  3  0  C 

Enter  the  numDer  cl  target  nuclei  for  the  113ln  sample: 
SE2  0 

Enter  total  oDserved  and  bkgnc  counts  in  the  392  keV  peak 


2404  C 

Enter  C,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 

i. 

Enter  tne  sell -absorption : 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  115ln  sample: 

3  0  C  6  C  3  0  C 

Enter  tne  number  of  target  nuclei  for  the  115lr:  sample: 
1E2 1 

Enter  total  observed  and  bkgnd  counts  in  the  336  keV  peak 
1  4  7  4  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  self-absorption: 

Enter  Tlrr,  TTrav  ano  TCnt  [sec]  for  the  135Ba  sample: 

300  60  300 

Enter  trie  number  of  taraet  nuclei  for  the  135Ba  sample: 
1E2  2 

Enter  total  observed  and  bkgnd  counts  in  the  268  keV  peak 
3357  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
i 

Enter  the  se 1 f -absor pt i on : 


Enter  TIrr,  TTrav  and  TCnt  [sec]  tor  the  179Hf  sample: 

3  0  C  6  C  3  0  C 

Enter  tne  number  of  target  nuclei  tor  the  179H£  sample: 
1E2C 

Enter  total  observed  and  bkqnd  counts  in  the  214  keV  peak: 
175  2  0 

Enter  C,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 

“I 

Enter  the  se 1 i -absorpt i on : 

Enter  TIrr,  TTrav  and  TCnt  [sec]  for  the  195Pt  sample: 

3  0  C  6  0  30L 

Enter  the  number  of  target  nuclei  for  the  195Pt  sample: 

Enter  total  observed  and  bkqnd  counts  in  the  130  keV  peak: 
109  9  0 

Enter  0,  1  or  1  for  a  Nal  solid,  well,  or  H'PGe  detector: 
Enter  the  se 1 f -absorpt i on : 


The  values  of  gateway  energy,  flux  and  standard  deviaton  are: 


2 . 6  6  0  0  0  C 

1 . 551038E+09 

2 . 143701E+07 

2 .790000 

1 . 398018E+09 

1.15740  5E  +  0  7 

1.980000 

2 . 621893E+09 

2 . 803883E+07 

2 .830000 

1 . 376944E+09 

2 . 34  52  8  5E  +  07 

1 .750000 

3 . 031389E+09 

0 . 000000E  +  00 

4.120000 

3 . 349735E+08 

5630891 . 000000 

1 .620000 

3 . 255730E+09 

0 . 000000E+00 

4.260000 

2 . 724333E+  06 

5 . 053730E+07 

1.960000 

2 . 772004E+0S 

5 . 653620E+07 

2 .460000 

1 . 743112E+09 

4 . 540216E+07 

2 .950000 

1 . 267657E+09 

0 . OOOOOOE+OO 

4 . 210000 

2 . 51447 5E+08 

6956013 . 000000 

2 .080000 

2 . 526803E+09 

6 . 036765E+07 

1 . 820000 

3 . 09  C  24  5E  +  09 

9 . 321678E+07 

C : \NEWFORT\ PROGRAMS  > 
C : \NEWFORT\ PROGRAMS  > 
C : \NEWF0RT\ PROGRAMS  > 
C : \NEWF0RT\ PROGRAMS > 
C: \ NEWPORT \PROGR AMS > 
C : \NEWFORT\PROGRAMS > 
C : \NEWFORT \ PROGRAMS  > 
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EXAMPLE 


C  Me  V 


DATA  GENERATION  PROGRAM  -  (ENFOLD)  -  Produces  observed  counts  represent 
mo  activations  present  m  hypothetical  samples  due  to  a  particuia 
given  input  x-ray  spectrum 

General  Instructions 

Gnomes  cf  input  spectrum,  detector  efficiencies,  type  of  random 
deviates,  and  spectral  plot  (if  desired)  are  made  through  menus. 

Spectrum,  menu  -  L 1  nac  ,  Ar  b  i  tr  a  r  y  or  Aurora  -  Alt-S 
Detector  menu  -  Nsl  veil  or  solid  or  HPGe  well  -  Alt-D 
Rand  on  deviates  menu  -  Exact  or  Poissor.  -  Alt-R 
Spectral  plot  menu  -  Linear, Log  or  MeV/MeV  -  Alt-P 


Spectrum :  Li nac  EEnd[Me”j=  6 


Detector:  Isomer  1: 

3"  x 

*5  !» 

Nal 

well 

I  some r  2 : 

2 "  >: 

?  M 

Na  I 

well 

I s omer  2  : 

3"  >: 

Nal 

well 

I s ome  r 4  ; 

2"  >: 

“>  1* 

Nal 

well 

Rant  dev:  Poisson 


Spectral  Distribution 


Energy [ MeV ] 

Phi [ l/cm~2-keV 

C  .  125 

3 . 68E+07 

0 .375 

1 . 20E+09 

0.625 

4  .  22E  +  09 

0.675 

4  . 18E  +  09 

1.125 

4  .  11E  +  09 

1.375 

3 . 88E+09 

1.625 

3  .  45E  +  09 

1.375 

2 . 65E+  09 

2 .125 

2  . 17E  +  09 

2 . 375 

1 . 81E+09 

2.625 

1 . 52E+09 

2.875 

1 . 29E  +  09 

3.125 

1 .10E  +  09 

3 . 375 

9 . 49E+08 

3 .625 

8  .  21E  +  08 

3.875 

7 . 14E+08 

4  .125 

6 . 23E+08 

4.375 

5 . 25E  +  08 

4 .625 

4 . 42E+08 

4.675 

3 . 55E+08 

5 . 125 

2 . 73E+08 

5.375 

2 . 01E+08 

5.625 

1 . 41E  +  08 

5.875 

9 . 42E+07 

Target  Isomer  Information 
Single  gateway  nuclei 


Isomerl-  87Sr  Isomer2-  89Y 


NTargetr 

=  1.0E+2G 

Ntargets 

=  1.0E+22 

EGare i heV ]  SG [ 

1C ' -29cm~2-keV] 

EGatetMeV]  SG  [  10~-29crr.  ~ 

i.  .  6  8 

2  0  06 

2 . 79 

592 

0 

c 

0 

0 

Q 

0 

r, 

u 

0 

0 

Ci 

0 

0 

C 

0 

0 

G 

0 

0 

0 

ExactAct [ 1/sec  ]  = 

2 . 98E-17 

ExactAct[l/sec]= 

8 . 06E-18 

Tl/2 [ s  ]  = 

1 . 01E+04 

Tl/2 [ s  ]  = 

1 . 61E  +  01 

IrrTls)  = 

200 

IrrTIs]  = 

300 

TTls] 

6  0 

TT  [  s  ] 

60 

CntTts]  = 

300 

CntTts]  = 

300 

DecCor r  i  = 

2 . 03E-02 

DecCorrl= 

7 . 55E-02 

DecCor r  2- 

9 . 90E-01 

DecCorr  2  = 

7 . 74E-02 

DeadTCor  = 

1 . 00E+0C 

DeadTCor  = 

1 . 00E+00 

E 1  [  k  e  V  ]  = 

388.4 

EltkeV]  = 

909 .15 

Intens  - 

0.82  3 

Intens  - 

0 .9914 

SelfAbs  - 

1.000 

SelfAbs  = 

1 . 000 

DetEff  = 

0  .339 

DetEff  = 

0 . 102 

ObsCtsl  = 

5144 

ObsCtsl  = 

14320 

StandDev= 

7  2 

StandDev= 

120 

E2 [ keV ]  = 

0 

E  2  [  k  e  V  ]  = 

0 

Intns  = 

0 

Intens  = 

0 

SelfAbs  = 

0 . 000 

SelfAbs  = 

0 . 000 

DetEff  = 

0.000 

DetEff  = 

0 . 000 

ObsCts2  = 

0 

ObsCts2  = 

0 

StandDe v= 

0 

StandDev= 

0 
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Isomer  3-  13  7Ba 

NTaraets=  I.OE^IS 


EGate [ Me” ] 

SG  [  10  '  -  29cn."  2-keV  ] 

1.9  c 

220  3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Exart Art [  1 

/  s  e  c  i  -  r.3GE-17 

Tl/I!s,  = 

1 . 53E+02 

I  r  r  T  i  s  J  = 

300 

TT  [  s  ] 

60 

CntT [ s  ]  = 

300 

DecCorr 1- 

5 . 66E-01 

DecCor  r  2  = 

5 . 47E-01 

DeadTCor = 

1 . 00E+00 

EICkeV]  = 

661  .66 

Intens  = 

0.901 

SelfAbs  = 

1 . 000 

DetEff  = 

0.181 

ObsCtsl  = 

8027 

StandDev= 

9  0 

E  2  t  k  e  V  i  = 

0 

Intens  - 

0 

SelfAbs  = 

0.000 

DetEff  = 

0.000 

ObsCtsl  = 

0 

StandDev= 

0 

Isomer4-  199Hg 
Ntargets=  1.0E+19 

EGate  [  Me  V  ]  SG  [  1 0 '  -- 29cm ~  2 -ke V  ] 


2.83 

3698 

0 

0 

0 

0 

(i 

0 

0 

0 

0 

0 

ExactAct [ 1/sec  ]  = 

4 . 91E-17 

Tl/2 [ s  ]  = 

2 . 5  6E  +  0  3 

IrrTls]  = 

300 

TTls] 

60 

CntTts)  = 

30  0 

DecCorr 1= 

7 . 69E-02 

DecCorr  2  = 

9 . 60E-01 

DeadTCor = 

1 . 00E+00 

EltkeV]  = 

158 . 376 

Intens  = 

0.53 

SelfAbs  = 

1 . 000 

DetEff  = 

0.577 

ObsCtsl  = 

3390 

StandDev= 

58 

E2 [ keV  ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0 . 000 

DetEff  = 

0 . 000 

ObsCts2  = 

0 

StandDev= 

0 
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Target  Isomer  Information 
Multiple  gateway  nuclei 


Isomeri-  77Se 
NTar gets=  1.0E+21 


isomeri-  lllCd 
Ntargets=  1.0E+20 


EGatelMe”;  SG [ 10 29cm ~ 2 -xeV ]  EGatelMeV]  SG[ 10~-29cm~2-keV j 


*1  • 


7  5 
12 

0 


0 

0 


2e9 

2891 

0 

0 

0 

0 


1.62 
4  .26 
0 
0 
0 
0 


529 

0 

0 

0 

0 


ExactAct[l/sec]= 

2 . 66E-17 

Tl/2 [ s  i  - 

1 . 75E+01 

IrrT [ s ]  = 

300 

TT  [  s  ] 

60 

CntTtsj  ^ 

300 

DecCor r  1  = 

9 . 22E-02 

DecCorr  2  = 

8 . 39E-02 

DeadTCor  = 

1 . OOE+OO 

EltkeV]  = 

161.92 

Intens 

0.524 

Self Ads  = 

1 . 000 

DetEff  = 

0 . 572 

ObsCtsl  = 

18502 

StandDe v= 

136 

E  2  [  k  e  V  3  = 

0 

Intns  = 

0 

SelfAbs  = 

0 .000 

DetEff  = 

0 .000 

ObsCts  2  = 

0 

StandDev= 

0 

ExactAct[l/sec]= 

4 . 79E-18 

Tl/2 [ s  ]  = 

2 . 9  2E  +  0  3 

I rr T [ s  ]  = 

300 

TT  [  s  ] 

60 

CntTIs]  = 

300 

DecCorr 1= 

6 . 79E-02 

DecCorr 2= 

9 . 65E-01 

DeadTCor = 

1 .OOE+OO 

EltkeV]  = 

245.384 

Intens  = 

C  .  94 

SelfAbs  = 

1.000 

DetEff  = 

0.472 

ObsCtsl  = 

4097 

StandDev= 

64 

E2 [ keV ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0 . 000 

DetEff  = 

0 .000 

ObsCts2  = 

0 

StandDev= 

0 
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Isomer  5-  1131  n 


NTarq 

ets=  5.0E+2G 

EGate [ Me V ] 

SG [ 1 0  ^  —  29cm' 

1  .96 

80 

4 . 51 

79" 

0 

0 

0 

Ci 

c 

0 

0 

0 

ExactAct [1/sec 

)=  5.81E-16 

Tl/2 [ s  3  = 

5 . 97E+03 

IrrTts]  = 

300 

TTis] 

60 

CntTts]  = 

300 

DecCor r 1= 

3 . 40E-02 

DecCor r  2- 

9 . 83E-01 

DeadTCor  = 

1 . OOE+OO 

El [ ke V  ]  = 

391.69 

Intens  = 

0.642 

SelfAbs  = 

1 . 0  0  c 

DetEff  = 

0  .337 

ODsCtsl  = 

6314 

StandDev= 

79 

E2 [ keV ]  = 

0 

Intens  = 

0 

SelfAbs  - 

0.000 

DetEff  = 

0.000 

ObsCts  2  = 

0 

StandDev= 

0 

lsomer4-  115In 
Ntarqets=  1.0E+21 

EGate  [  HeV  ]  SG  [  10~-29cnT  2-keV  ) 


2 .46 

123 

4 . 82 

1228 

0 

0 

0 

0 

0 

0 

0 

0 

Exact Act [ 1/sec  3  = 

6 . 69E-18 

Tl/2  1  s  ]  = 

1 . 57E+04 

IrrTls]  = 

300 

TT  [  s  ] 

60 

CntTts]  = 

300 

DecCorr 1= 

1 . 32E-02 

DecCor r  2- 

9 . 93E-01 

DeadTCor  = 

1 . 00E+00 

EllkeV]  = 

336.26 

Intens  = 

0 .458 

SelfAbs  = 

1 . 000 

DetEff  = 

0.383 

ObsCtsl  = 

4726 

StandDev= 

69 

E2 [ keV ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0 . 000 

DetEff  = 

0 .000 

ObsCts2  = 

0 

StandDev= 

0 
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Target  Isomer  Information 
Multiple  gateway  nuclei 


IsomerS-  135Bs 
N’Target?-  1.0E-r22 

EGa  t  e  t  Me  V  j  S G  [  1 0 '  -  2  9  cm '  2  -  k  e  V  ] 


O  Q  v 

211 

4.21 

2110 

0 

0 

0 

0 

0 

0 

0 

0 

ExaotAct [ 1/sec  ]  = 

1 . 52E-17 

Tl/2ls]  = 

1 . 03E+05 

1 rr T [ s  ]  = 

300 

TT  ( s  ] 

60 

CntTIs]  = 

300 

DecCorr 1= 

2 . 01E-03 

DecCorr  2  = 

9  .  99E-C1 

DeadTCor  = 

1 . 0GE+00 

EllkeV]  = 

266.27 

Intens  = 

0.156 

SelfAbs  = 

1.000 

BetEff  = 

0.448 

ObsCtsl  = 

6371 

StandDev= 

80 

E2 [ keV  J  - 

0 

Intns  = 

0 

SelfAbs  = 

0.000 

BetEff  = 

0 . 000 

ObsCts2  = 

0 

StandBev= 

0 

Isomer6-  179Hf 
Ntargets=  1.0E+20 

EGare [  Me  V  ]  SG [ 1 0 ~ - 29cm' 2 -ke V ] 


2 .06 

494 

4.69 

4942 

0 

0 

0 

0 

0 

0 

0 

0 

ExactAct[l/sec]= 

3.18E-17 

Tl/2 [ s ]  = 

1 . 87E+01 

IrrTl s ]  = 

300 

TT  [  s  ) 

60 

CntTIs]  = 

300 

DecCorr 1= 

1 . 08E-01 

DecCorr  2  = 

8 . 98E-02 

DeadTCor  = 

1 . OOE+OO 

El [ ke V ]  = 

214.31 

Intens  = 

0.952 

SelfAbs  = 

1.000 

DetEff  = 

0 . 507 

ObsCtsl  - 

4474 

StandDe v= 

67 

E2 [ keV ]  = 

0 

Intens  = 

0 

SelfAbs  = 

0 . 000 

BetEff  = 

0 .000 

ObsCts2  = 

0 

StandDev= 

0 
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Isomer"-  195Pt 
NTargets=  5.0E+22 

EGate [ MeV ]  SG [ 1 0 29cm" 2-keV ] 


1.82 
5 . 4  8 
0 
0 
0 
0 


0 

0 

0 

0 


ExactAct [ 1/sec  ]  =  1.01E-17 


Tl/2[s]  = 
IrrTEs)  = 
TT  E  s  ] 

CntT [si  = 
DecCorr  1- 
DecCor r  2  = 
DeadTCor = 

EllkeV]  = 
Intens  = 
SelfAbs  = 
DetEff  = 

ObsCtsi  - 
Stand be v= 


2 .  47E  +  05 
300 
60 
300 
5 . 98E-04 
I . 00E+00 
1 . 00E+00 

129.76 
0 .02885 
1.000 
0.616 


41 


E  2  [  k  e  V  i  =  0 

Intens  =  0 

Se 1 f Ads  =  0.000 

DetEff  -  0.000 

ObsCtsi  =  0 

Standbev=  0 
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C:  \N2WF0ET \PROGRAMS  .'unfold 
Enter  1  for  a  c w  exposure: 

1 

Enter  tnt  end  point  energy  [MeV]  of  the  x-ray  device: 

C 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  87Sr  sample: 

300  60  30C 

Enter  the  numoer  of  target  nuclei  for  the  87Sr  sample: 

IE  2  0 

Enter  total  observed  and  bkgnd  counts-  in  the  388  keV  peak: 
:14  4  0 

Enter  0,  1  or  1  for  a  Nal  solid,  veil,  or  HPGe  detector: 
Enter  the  self-absorption: 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  89Y  sample: 
30C60300 

Enter  the  numoer  of  target  nuclei  for  the  89Y  sample: 

1E2  2 

Enter  total  observed  anc  bkgnd  counts  in  the  909  keV  peak: 
143200 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  sel f -absorpt i on : 


Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  137Ba  sample: 

3  0  C  60  300 

Enter  the  number  of  target  nuclei  for  the  137Ba  sarple: 

1E19 

Enter  total  observed  and  bkqnd  counts  in  the  661  keV  peak: 

3027  0 

Enter  0,  1  or  1  for  a  Nal  solid,  veil,  or  HPGe  detector: 

Enter  the  se 1 f -absorpt i on : 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  199Hg  sample: 

300  60  300 

Enter  the  number  of  target  nuclei  for  the  199Hg  sample: 

1E15 

Enter  total  oDserved  and  bkqnd  counts  in  the  158  keV  peak: 

2  39  0  0 

Enter  0,  1  or  2  for  a  Nal  solid,  veil,  or  HPGe  detector: 

1 

Enter  the  se 1 f -absorpt i on : 

1 

Enter  1  to  examine  the  fit  employed  in  the  single  gateway  region 


The  coefficients  of  the  power  series  are: 
n  =  1,  A  =  2.942580E+09 

=  2,  A  =  1.75878  6E+  09 

n  =  3,  A  =  -1 . 497582E+09 

r,  =  4,  A  =  2 . 387671E  +  08 

Chi  squared  =  8.092223E-01 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  77Se  sample: 


300  6 C  300 

Enter  the  number  cf  tarqet  nuclei  for  the  77Se  sample: 

1E2  1 

Enter  total  observed  and  Lxqnd  counts  in  the  162  keV  Deal. 
185C2  0 

Enter  G,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  tne  self-absorption: 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  lllCd  sample: 

300  50  30C 

Enter  tne  number  of  target  nuclei  for  the  lllCd  sample: 
1E2C 

Enter  total  observed  and  bkqnd  counts  in  the  240  keV  peak 
409"  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 

J. 

Enter  the  se 1 f -abs or pt i on : 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  113In  sample: 

300  60  30C 

Enter  the  number  of  target  nuclei  for  the  1131n  sample: 
5E2  0 

Enter  total  observed  and  bkgnd  counts  in  the  392  keV  peak 
6314  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 

1 

Enter  tne  se 1 f -abs or pt i on : 

1 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  115ln  sample: 

3  0  C  6  C  300 

Enter  the  number  ox  target  nuclei  for  the  115In  sample: 

IE  2  1 

Enter  total  observed  and  bkgnd  counts  in  the  336  keV  peak 
4726  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 

i. 

Enter  the  self-absorption: 

7 

Enter  Tlrr,  TTrav  and  TCnt  [sec]  for  the  135Ba  sample: 

300  60  300 

Enter  the  number  of  target  nuclei  for  tne  135Ba  sample: 
1E2  2 

Enter  total  observed  and  bkgnd  counts  in  the  260  keV  peak 
6371  0 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 


Enter  the  se 1 f -absor pt i on : 


•  *" 
i 

Enter  TIrr,  TTrav  and  TCnt  [sec]  for  the  179Hf  sample: 

3  0  C  6  C  3  0  C 

Enter  the  number  of  target  nuclei  for  the  179Hf  sample: 

1E2 1 

Enter  total  observed  and  bkgnd  counts  in  tne  214  keV  peak: 

4  4  *7  4  1 

Enter  G,  1  or  2  tor  a  Nal  solid,  well,  or  HPGe  detector: 

i 

Enter  the  self -absorption: 

Enter  TIrr,  TTrav  and  TCnt  [sec]  for  the  195Pt  sample: 

3  0  C  6  C  3  0  0 

Enter  the  number  of  target  nuclei  for  the  195Pt  sample: 

5E2  2 

Enter  total  observed  and  bkgnd  counts  in  the  130  keV  peak: 
16  5  2  C 

Enter  0,  1  or  2  for  a  Nal  solid,  well,  or  HPGe  detector: 
Enter  the  se 1 f -absor pt i on : 


The  values  of  gateway  energy,  flux  and  standard  deviaton  are: 


2 .660000 

1 . 524077E+09 

2 . 12498 7E+07 

2 . 790000 

1 . 3  7  214  7E+  09 

1 . 14664  5E  +  07 

1.980000 

2 . 406901E+09 

2 . 686467E+07 

2.830000 

1 . 35417 5E+09 

2 . 325813E+07 

1 .750000 

2 . 713751E+09 

0 . OOOOOOE+OO 

4.120000 

6 . 482090E+08 

6759861 .000000 

1.620000 

2 . 876682E+09 

0  .  OOOOOOE  +  OO 

4.260000 

6 . 003942E+08 

4 . 566513E+07 

1.960000 

2 . 4  3  4  4  9  3E  +  09 

0 . 000000E+00 

4.510000 

4 . 8  6  4  2  9  4E+  0  8 

1. 529281E+07 

2 . 460000 

1 . 7  6  09  3  3E  +  0  9 

0  .  OOOOOOE  +  OO 

4.820000 

3 . 834140E+06 

1 . 77  51 6  6E+  07 

2 .950000 

1 . 228004E+09 

0 . OOOOOOE+OO 

4.210000 

5 . 949821E+08 

9289119 .000000 

2.080000 

2 . 270358E+09 

0  .  OOOOOOE  +  OO 

4.690000 

4 . 180527E+08 

1 . 031957E  +  07 

1 . 820000 

2 . 622402E+09 

0 . OOOOOOE+OO 

5 . 480000 

2 . 020995E+08 

1 . 907023E+07 

C : \NEWFORT\PROGRAMS  > 
C : \NEWFORT\ PROGRAMS  > 
C  :  \NEWFORT\ PROGRAMS > 
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